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ABSTRACT 



In recent years significant attention has been attracted to proposals which 
utihze DNA for nanotechnological apphcations. Potential applications of these 
ideas range from the programmable self-assembly of colloidal crystals, to biosensors 
and nanoparticle based drug delivery platforms. In Chapter I we introduce the 
system, which generically consists of colloidal particles functionalized with specially 
designed DNA markers. The sequence of bases on the DNA markers determines 
the particle type. Due to the hybridization between complementary single-stranded 
DNA, specific, type-dependent interactions can be introduced between particles 
by choosing the appropriate DNA marker sequences. In Chapter 11 we develop 
a statistical mechanical description of the aggregation and melting behavior of 
particles with DNA-mediated interactions. A quantitative comparison between the 
theory and experiments is made by calculating the experimentally observed melting 
profile. In Chapter III a model is proposed to describe the dynamical departure 
and diffusion of particles which form reversible key-lock connections. The model 
predicts a crossover from localized to diffusive behavior. The random walk statistics 
for the particles' in plane diffusion is discussed. The lateral motion is analogous to 
dispersive transport in disordered semiconductors, ranging from standard diffusion 
with a renormalized diffusion coefficient to anomalous, subdiffusive behavior. In 
Chapter IV we propose a method to self-assemble nanoparticle clusters using 
DNA scaffolds. An optimal concentration ratio is determined for the experimental 
implementation of our self-assembly proposal. A natural extension is discussed 

XV 



in Chapter V, the programmable self-assembly of nanoparticle clusters where the 
desired cluster geometry is encoded using DNA-mediated interactions. We determine 
the probability that the system self-assembles the desired cluster geometry, and 
discuss the connections to jamming in granular and colloidal systems. In Chapter 
VI we consider a nanoparticle based drug delivery platform for targeted, cell specific 
chemotherapy. A key-lock model is proposed to describe the results of in-vitro 
experiments, and the situation in-vivo is discussed. The cooperative binding, and 
hence the specificity to cancerous cells, is kinetically limited. The implications for 
optimizing the design of nanoparticle based drug delivery platforms is discussed. 
In Chapter VII we present prospects for future research: the connection between 
DNA-mediated colloidal crystallization and jamming, and the inverse problem in 
self-assembly. 
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CHAPTER I 



A DNA-COLLOIDAL PRIMER 
1.1 Miniaturization 

Advances in science have made possible the manipulation of matter on a smaller 
and smaller scale. Controlling the spatial arrangement of atoms and molecules 
enables the control of bulk material properties. Miniaturization has attracted 
significant attention in its own right, particularly with respect to integrated circuit 
design for computer hardware. This trend, known as Moore's Law, states that 
the number of transistors which can be placed on an integrated circuit has been 
increasing exponentially, approximately doubling every two years [5]. Independent 
of our ambition to quench the thirst for increased computing power, miniaturization 
will likely play an important role in the future of medical science. The ability 
to engineer nanodevices which interact with individual cellular components has 
a number of potentially exciting applications, ranging from smart drug delivery 
vehicles [6], [4], [7], [8], [9], [10] to biosensors which can detect an astonishingly 
low concentration of pathogens [11], [12], [13]. The realization of these goals 
depends fundamentally on our ability to control the structure and arrangement of 
individual components on the nanoscale. On these lengths we encounter problems 
with traditional top-down assembly approaches to miniaturization, for example 
lithography [14]. One proposed resolution is to proceed from the bottom- up, 
harnessing the incredible molecular recognition properties of DNA [15], [16], [17], 
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[18]. 

1.2 DNA 

Deoxyribonucleic acid, hereafter simply DNA, is a biopolymer which contains 
the genetic information for the function of all living organisms. The macromolecule 
consists of a sugar-phosphate backbone chain with the saccharide unit carrying a 
nucleotide of four possible types [19]. The primary structure of DNA refers to 
the sequence of these nucleotides from the four letter DNA alphabet consisting 
of cj^osine (C), thymine (T), adenine (A), and guanine (G). The secondary 
structure of DNA refers to the short range order which manifests itself as a result 
of interactions between monomers which are in close proximity [20]. Hydrogen 
bonding between complementary DNA base pairs results in a DNA double hehx, 
in which two DNA molecules wind around each other. The complementarity rule 
states that adenine bonds with thymine, and cytosine bonds with guanine. The 
double helix is approximately 2nm in diameter, and the repeat in the direction of 
the helix axis is every 3AA which is about every 10 base pairs. The energy gain 
associated with forming a base pair in the double helix is comparable to the hydrogen 
bond energy AE ~ O.leV. The formation of the T-A (C-G) pair is a result of two 
(three) hydrogen bonds. As a result the characteristic energy required for double 
stranded DNA to denature and form two single strands is comparable to the thermal 
energy at room temperature T ~ 300K. In many respects DNA appears to be an 
excellent candidate to control matter on the nanoscale. The interactions between 
nucleotides are highly specific. In addition, the number of potential sequences 
4^ = exp(A^log4) grows exponentially with the number of nucleotides N. 
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1.3 Polymer Physics 

Part of the usefulness of DNA in controlling matter on the nanoscale stems 
not from its chemical specifics, but general conformational properties of long 
chain-hke molecules [21]. Here we introduce some of the basic ideas in studying the 
conformations of polymers which will be of use later on. 

A very idealized model of a polymer is a sequence of N rigid links of length /, 
where the direction between consecutive links is independent. In this freely-jointed 
model the end to end distance of the polymer chain r can be expressed in terms of 
the bond vectors Uj = Xj+i — Xj where Xj is the radius vector of the i^'^ segment. 

N 

r=J]u, (1.1) 

i=l 

The radius of gyration of the chain Rg is defined in terms of the average mean-squared 
displacement. 

/ N N \ N 

= (r • r) = ( E E ) = E <"^'> + E ("^•^^■) = (1-2) 

\ i=l j=l I i=l i^j 

The cross terms vanish when averaged since we assumed the angular orientation of 
the links was uncorrelated. Note that the characteristic size of the ideal polymer 
Rg — N^l is significantly smaller than that of the fully extended chain. In the limit 
of large N the probabihty distribution function PAr(r) for a particular end to end 
distance r is Gaussian. 

This statement follows from the central limit theorem, since the end to end distance 
can be expressed as a sum of independent bond vectors. Alternatively one can 
consider the polymer configuration as a random walk, in which case -Piv(r) satisfies 
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the diffusion equation [20], [22]. Hence at fixed r the entropy of the polymer chain is 

^(r) = log(P^(r)) = const - ^ (1.4) 

Note that thourhgout this thesis we refrain from writing the Boltzmann constant, 
choosing natural units with /c^ = 1. Since there is no interaction energy in this 
model the free energy can be written as follows: 

3T 

F(r) ^E- TS{v) = const + — -r^ (1.5) 



If we stretch the chain by applying a stretching force f on both ends of the 
polymer the free energy increases. In equilibrium the corresponding elastic restoring 
force fei — — f . The extension of the polymer chain R as a result of applying a 
stretching force is 

ST 

f = VF(r)U^=^R (1.6) 

which is vahd provided the chain is not stretched too much |R| <^ Nl. Hence we 
see that the ideal polymer behaves like a mechanical spring with spring constant 
k — 3T/ {Nl'^). The chain stretches along the direction of the apphed force, and the 
corresponding restoring force is of purely entropic origin (i.e. since there are fewer 
configurations of the stretched chain). 

On long enough length scales the single chain will be ideal. When excluded 
volume interactions are included between the monomers, we expect to see deviations 
from the ideal chain behavior. Flory presented the following argument [23] to 
determine how the size of the polymer chain depends on the number of monomers 
N. We expect that the excluded volume between monomers will favor swelling of 
the chain. If chain is confined to a volume the average monomer concentration 
c ~ N/ B?. As a result the total repulsive energy associated with monomer-monomer 
interactions is proportional to Fj-^p ~ Tv(?B? where we have introduced the excluded 
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volume parameter v which in general may be temperature dependent. However, 
stretching the chain costs entropy, so there is a contribution to the free energy 
Fei ~ Minimizing the total free energy F — Frep + F^i to determine the 

preferred chain size Rp we have 

^ (F +F) - ^Z:!^,2Ti?_ 

Rf ~ Nhiv^^Nh (1.8) 

The Flory exponent = 3/5 gives the dependence of the chain size Rp '-^ on the 
monomer number N . In writing the last equality I have estimated the excluded 
volume parameter v ~ Note that the chain is stretched as compared to the ideal 
chain which has v — 1/2. 

With this in mind, we can return to the question of determining the stretching 
response of a chain with excluded volume interactions. Here we present a scaling 
argument due to Pincus [24]. The characteristic length which enters the problem is 
the flory radius Rp — N'^l- The other parameters of the problem are the magnitude 
of the stretching force / and the thermal energy T. A scahng function (p{x) with 
dimensionless argument x = {RFf)/T is introduced to determine the elongation of 
the polymer R in response to the stretching force. 

R ~ Rpifiix) (1.9) 

When the stretching force is small the chain is weakly perturbed and the response 
should be proportional to /. Hence for a; <C 1 we have (p{x) ~ x and 

R Rfx ^ ml (^-pj . (1.10) 

So we see that the spring constant of the chain with excluded volume interaction 
k ~ N~^/^ is reduced as compared to that of the ideal chain k ~ N~^. In the 
opposite regime of strong stretching we require that the extension R be linear in N . 
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Hence for x ^ 1 we assume (p{x) ~ x"^ and determine m — 1 — {1/ u) = 2/3 which 
satisfies this condition. 



For strong stretching, the chain with excluded volume interaction has a nonlinear 
force-extension relation / ~ i?^/^ which deviates from the linear Hooke's Law for 
the ideal polymer chain. These results will be of interest in later chapters when we 
model the interaction of colloids which are connected by polymer springs. 

In a real polymer system there will be correlations between adjacent links, in 
which case our assumption (uj-Uj) = in Eq. 1.2 is no longer valid. The persistence 
length Ip of the polymer chain provides a measure of the chain fiexibility, and is 
roughly the maximum length for which the polymer chain remains straight. Let 
9{s) be the angle between two segments of the chain separated by a distance s. In 
these terms the persistence length is defined as [20] 



The persistence length of single-stranded DNA [Ip ~ Inm) is significantly shorter 
than that of double-stranded DNA {Ip ~ 50nm). The double helix structure is quite 
rigid, whereas the single strand is more flexible. For lengths L < Ip the chain can 
be treated effectively as a rigid rod. 



Here we present one approach whereby DNA can be used to organize particles 
on the nanoscale. The general idea is to graft many DNA strands onto the surface 
of a colloidal particle [25] . The size and chemical composition of the colloid depends 
on the application. In some experiments polystyrene beads with diameter d ~ l/xm 
are utilized for this purpose [1], [26]. Another common experimental approach [27], 
[28], [29] is to use gold nanoparticles with d ~ lOnm. In this case the grafting is 
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(1.11) 




(1.12) 



1.4 DNA Grafted Colloids 
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made possible by attaching a thiol group to one end of the DNA strand which binds 
to the surface of the gold nanoparticle. The result is a system of monodisperse 
"octopus-hke" particles where each particle has many DNA arms. One end of each 
DNA chain is attached to the surface of the particle, and the other end is free. 
In preparing such a system the experimenter can control both the average DNA 
grafting density, and the particular nucleotide sequence of the DNA arms. Note 
that preparing these "octopus-hke" particles relies on diffusion of DNA chains which 
adsorb to the particle surface. This adsorption process is random or stochastic, as 
a result one cannot control the exact number N of DNA chains attached to a given 
particle. Instead one controls the average number of DNA chains per particle (N) 
by choosing the appropriate ratio of the total DNA concentration Cdna to the total 
particle concentration C particle during preparation. 

{N) = (1.13) 

^particle 

This parameter {N) completely defines the probability distribution for the number 
of DNA arms N attached to a given particle, which due to the random character of 
the preparation process must have the Poisson form. 

Here P{N) is the probability that a particle has exactly N DNA arms, with (N) the 
average number of DNA arms on the particles. 

The ability to control the sequence of DNA nucleotides attached to the particles 
leads to interactions between particles of different types in solution. We say 
that the "type" or "color" of the particle is determined by the sequence of DNA 
nucleotides attached to the particle. For example, consider particles grafted with 
many single-stranded DNA with sequence ACTGAG. Wc call these "red" particles. 
We could also prepare "green" particles with sequence CTCAGT. Here I label the 
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sequences with the following rule. The first letter is the base which is closest to 
the grafting point, out to the last letter which is the base at the free end of the 
DNA chain. Note that I have chosen the green sequence complementary to the 
red one as dictated by the rule for complementary hydrogen bonding. In solution 
when DNA arms of the red particles encounter DNA arms of the green particles 
these two single strands of DNA can hybridize to form a double strand. Provided 
that we are working under appropriate experimental conditions (temperature, salt 
concentration, etc.) the formation of the double strand results in a lower free energy 
state than if the two strands were denatured. This provides a practical method 
to hnk particles through the formation of a DNA bridge. The bond that results 
between particles connected by DNA bridges is reversible, since we can change the 
temperature or pH of the solution to denature the two DNA strands composing the 
bridge. The binding energy for the formation of a DNA bridge will depend on a 
number of factors, including the length of the complementary DNA sequence and 
properties of the DNA chains attached to the particles [30] . For now we simply note 
that the interaction is highly-specific and tunable. 

1.5 Interactions 

The DNA-colloidal system we are considering is quite complex. In general 
the interaction potential between colloids in solution combines specific (or type- 
dependent) interactions with non-specific (type-independent) interactions. The 
specific interactions pertain to the formation of DNA bridges between colloids as a 
result of DNA-DNA hybridization. The specificity is determined by the sequence 
of DNA nucleotides attached to the particles and the complementary rule for DNA 
base pairing. The non-specific interactions include all the interactions which are 
independent of the particular DNA sequence. For example, this includes the van der 
Waals attraction and electrostatic repulsion as described by the DLVO theory [31], 



8 



[32]. In addition we must take into account the steric repulsion between colloids that 
arises from grafting many DNA chains to the surface of the particles. At first glance 
a quantitative treatment of the system appears discouraging given the complexity 
and diversity of the interactions. However, by comparing the characteristic energy 
and length scales we will see that the most important interactions for our purposes 
are those directly related to DNA, specifically DNA-DNA hybridization and steric 
repulsion. 

We first consider the non-specific interactions between colloids in solution. The 
electrostatic interactions between charged colloids in ionic solution are described by 
the Poisson-Boltzmann equation. Because the equation is nonlinear an analytic 
treatment is generally only possible with simple geometries in the context of some 
approximation scheme. In the context of the Debye-Hiickle approximation the 
equation can be linearized to obtain [33], [34], [35] the pair potential between two 
spherical colloids of radius a carrying fixed charge —Ze. 



Here Ib = e^/(eT) is the Bjerrum length and e is the dielectric constant of 
water [20]. For water at room temperature e ~ 80 which gives Ib ~ O.Jnm. 
The presence of counterions in solution leads to screening of the electrostatic 
potential. For monovalent counterions of concentration n the Debye screening 
length — 1/ y/i-TxlBn. The Debye length is the length at which the counterions 
screen out electric fields. For example, in a NaCl solution with concentration 
0.2M the Debye length k'^ k. 0.68 nm. This ion concentration is typical of many 
animal fluids. This estimate indicates that stabilizing colloidal suspensions against 
non-specific aggregation electrostatically is not a particularly appealing method due 
to the incredibly short range of the resulting repulsive potential. This is especially 
true in many biological applications where temperature and ion concentration are 
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not set by the experimenter. In general we will consider situations where the colloids 
themselves are not charged, and electrostatic interactions can be neglected. 

Even if the colloids are not charged, we still need to consider the DNA. In 
solution the phosphates which constitute the DNA backbone dissociate and each 
carries a negative charge. Because each of the links carries charge, we might 
expect that the repulsive interactions will lead to highly stretched conformations 
of the chain Rg ~ N. Here N is the number of monomers in a single DNA chain. 
However, in ionic solution the charges are screened. In fact for a strongly charged 
polyelectrolyte in ionic solution the counter ions condense on the chain, effectively 
neutralizing its charge [36]. Roughly speaking, the counterions condense once 
the linear charge density of the chain p exceeds the critical value p^rit = 
Electrostatic effects play a role in determining certain properties of the DNA chains, 
for example they increase the persistence length as compared to a neutral chain. 
However, from our perspective the fact that the DNA backbone is charged will not 
be of great importance. 

We now consider the van der Waals interaction between colloids. Consider an 
atom which on average has a spherically symmetric charge distribution. Quantum 
mechanical fluctuations of the valence charge give rise to an instantaneous dipole 
moment. The instantaneous dipole results in an electric field at a distance r from 
the atom E ~ This field induces a dipole moment "p^ ~ in a nearby 

atom. The resulting interaction energy U ~ —~p ■ E ~ — Assuming that 

the interaction between a collection of atoms is pairwise additive and nonretarded 
one can write [31] the following expression (in the Derjaguin approximation) for 
the interaction potential between two spheres of radius a. Here the spheres are 
separated by a surface to surface distance D — r — 2a and the expression is valid for 
D -C a. 



Uvdw{.r) 
T 



A a 



(1.16) 
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Here A — A/T is the reduced Hamaker constant. At T = 300K the reduced 
Hamaker constant yl ~ 3 for polystyrene in water, and ^4 ~ 76 for gold in water. 
For quantitative comparisons Eq. 1.16 is not particularly useful. A more detailed 
treatment is required which takes into account the effects of retardation. 

The resulting attraction is insignificant when compared to the specific attraction 
generated by DNA hybridization [37]. For example, in a recent study with 
micron sized polystyrene spheres, the van der Waals attraction was estimated to be 
Uvdw — — 3T at surface to surface separations of D = 14nm and U^dw — — lOT at 
D — lOnm. This is to be compared with the energy scale for the DNA hybdriziation, 
which will depend on the length of the complementary hybridization sequence. For 
a 15 base pair linker at room temperature Udna — — 30r per DNA bridge! 

We now consider the steric repulsion of the DNA chains which prevents the 
non-specific aggregation of colloids. Understanding the behavior of polymer brushes 
is an active field of research. Treatments of increasing complexity are available, 
from scaling arguments to self-consistent field theories [38], [39], [40], [41]. Here we 
present a simple argument to outline the qualitative behavior of grafted polymer 
brushes. As the surface grafting density a of the DNA chains increases, there is 
a competition between entropic and excluded volume effects. There is an energy 
penalty associated with monomer-monomer contacts which favors stretching of the 
chain. However stretching the chain costs entropy as discussed earlier. The result is 
the formation of a DNA brush on the surface of the colloid. These brushes interact 
giving rise to a repulsive potential between particles grafted with polymer chains. 

Writing the competition between the excluded volume and entropic interactions 
the free energy per chain F in the brush of height h is 



Here {Na)/h is the average monomer concentration in the brush with surface 




(1.17) 
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grafting density a and v is the excluded volume parameter. Minimization with 
respect to h gives the free energy per chain and the equilibrium brush height 
where we have estimated i> ~ Z^. 

~ TN{fa)^ (1.18) 
K ^ Nlifa)^^ (1.19) 

The resulting DNA brush is characterized by highly extended conformations of the 
DNA chains, in particular the equilibrium height of the brush is proportional to 
the number of monomers in a chain N. 

There is an energy penalty associated with compressing the brushes once the 
particle separation D < 2h^. A more detailed treatment of the problem takes into 
account the distribution of chain ends within the brush. By making the analogy to 
an associated quantum mechanical problem the authors of [40] have calculated the 
free energy penalty associated with compressing the brush to a height h < h^. Here 
we quote the result for the free energy per chain F{u) associated with compressing 
the brush to a height h < h^. The dimensionless parameter u — h/h^. 

(1.20) 

An order of magnitude estimate [29] for compressing the DNA brush below its 
equilibrium height repulsive gives several T per DNA chain. Therefore by tuning 
the brush height steric repulsion prevents particles from ever approach at separations 
close enough to feel a significant effect of the van der Waals attraction. Grafting 
polymers to the particle surface is a controlled technique one can utilize to prevent 
non-specific aggregation of particles in solution. Note that during this discussion we 
considered a DNA brush, but the mechanism is largely independent of the particular 
monomer chemistry. Another water soluble polymer could play a similar role, one 
common choice in experiments is polyethylene glycol (PEG). 



F{u) 



5F, 



1 9 M 

U 
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The result of this discussion indicates that in the DNA-coUoidal system we will 
consider, the pertinent interactions are those directly relating to the DNA. There 
is a specific attraction associated with DNA hybridization, and a non-specific steric 
repulsion which arises as a result of grafting many DNA chains on the surface of the 
colloids. 

1.6 Literature Review 

In this section we will highlight some of the literature which addresses problems 
related to the topics of this thesis. In the past two decades, there have been 
a number of experimental advances in DNA based self-assembly. These ideas 
originally stem from work in the lab of Ned Seeman, who introduced the first 
schemes for building nanostructured objects using specifically designed DNA [15], 
[18]. This approach has been adapted to demonstrate the ability of DNA to 
rationally assemble aggregates of colloidal particles. There have been and number of 
important contributions, including research in the groups of Mirkin [25] , Alivisatos 
[42], Soto [16], and many others [43], [44], [45], [46], [47], [48], [1], [26], [49], [50], 
[51], [11], [52], [53], [54], [37]. One particularly interesting recent advance is the 
self-assembly of colloidal crystals using DNA-mediated interactions by the groups 
of Gang [29] and Mirkin [28]. These systems have also attracted attention from a 
theoretical perspective. In one of the first theoretical works on the subject [55] , 
Tkachenko studied the equilibrium phase behavior for a binary system of particles 
decorated with DNA. The system exhibits a diverse spectrum of crystaUine phases, 
including the diamond phase which is of interest for the self-assembly of photonic 
crystals. 

Some previous theoretical work on the aggregation and melting behavior of 
DNA-coUoidal assemblies include the work of Jin et al. [2], Park and Stroud [27], 
and Lukatsky and Prenkel [56]. These authors studied the aggregation behavior 
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and optical properties of DNA-mediated colloidal assemblies. One drawback to 
the previous work is that the results were based on phenomenological or lattice 
based models which give limited insight into the physics underlying the aggregation 
phenomena. In Chapter II [30] we develop an off-lattice, statistical mechanical 
description of aggregation and melting in these systems. The results of the theory 
are compared quantitatively to recent experiments by the groups of Chaikin [1] and 
Crocker [26]. There are connections between this aggregation behavior and the 
sol-gel tranasition in branched polymers [21]. Other soft matter systems exhibit 
similar phenomena, for example a system of microemulsion droplets connected by 
telechelic polymers [57]. 

In addition to the work on bulk systems, DNA is a promising candidate to 
self-assemble small clusters of particles, or nanoblocks. Independent of the DNA 
based studies, Manoharan et al. [3] devised a scheme to self-assemble small clusters 
of microspheres. The microspheres are attached to the surface of liquid emulsion 
droplets, and the clusters self-assemble by removing fluid from the droplet. The 
clusters are packings of spheres that minimize the second moment of the mass 
distribution. This packing sequence is somewhat ubiquitous in soft matter systems. 
Glotzer et al. [58], [59] have demonstrated that cone-shaped clusters with N ^ 10 
particles self-assemble into clusters with the same packing sequence as [3]. This 
result is not necessarily expected, since the self-assembly processes are driven by 
different mechanisms. In the experiments capillary forces are responsible for the 
assembly process, whereas in the simulations the interactions between cone-particles 
are anisotropic and highly specific. Similar ideas can be used to explain the 
structure of prolate virus capsids [60]. In Chapter IV [61] we propose a method 
to self-assemble clusters of particles with the same packing sequence, where the 
self-assembly process is mediated by DNA. Other recent studies of the DNA based 
assembly of nanoscale building blocks include [62] and [63] . 
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CHAPTER II 



DNA-MEDIATED COLLOIDAL 
AGGREGATION 

2.1 Introduction 

In the past ten years, there have been a number of advances in experimental 
assembly of nanoparticles with DNA-mediated interactions [25], [64], [49], [50], 
[42], [51]. While this approach has a potential of generating highly organized and 
sophisticated structures [55], [65], most of the studies report random aggregation 
of colloidal particles [1], [26]. Despite these shortcomings, the aggregation and 
melting properties may provide important information for future development of 
DNA-based self-assembly techniques. These results also have more immediate 
implications. For instance, the observed sharp melting transition is of particular 
interest for biosensor applications [11]. For these reasons the development of a 
statistical mechanical description of these types of systems is of great importance. 
It should be noted that the previous models of aggregation in coUoidal-DNA systems 
were either phenomenological or oversimplified lattice models [2], [56], [27], which 
gave only limited insight into the physics of the phenomena. 

In this chapter [30] , we develop a theory of reversible aggregation and melting 
in coUoidal-DNA systems, starting from the known thermodynamic parameters of 
DNA (i.e. hybridization free energy AG), and geometric properties of DNA-particle 
complexes. The output of our theory is the relative abundance of the various 
colloidal structures formed (dimers, trimers, etc.) as a function of temperature, 
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as well as the temperature at which a transition to an infinite aggregate occurs. 
The theory provides a direct link between DNA microscopies and experimentally 
observed morphological and thermal properties of the system. It should be noted 
that the hybridization free energy AG depends not only on the DNA nucleotide 
sequence, but also on the salt concentration and the concentration of DNA linker 
strands tethered on the particle surface [66]. In this paper AG values refer to 
hybridization between DNA free in solution. 

In a generic experimental setup, particles are grafted with DNA linker sequences 
which determine the particle type(A or B). In this chapter we will restrict 
our attention to a binary system^. These linkers may be fiexible or rigid. A 
selective, attractive potential between particles of type A and B can then be 
turned on by joining the linkers to form a DNA bridge. This DNA bridge 
can be constructed directly if the particle linker sequences are chosen to have 
complementary ends. Alternatively, the DNA bridge can be constructed with the 
help of an additional hnker DNA. This additional linker is designed to have one 
end sequence complementary to the linker sequence of type A particles, and the 
other end complementary to type B. The properties of the DNA bridge formed will 
depend on the hybridization scheme(see figure 2.1). 

The plan for the chapter is as follows. In section 2.2 we provide a description of 
the problem. In section 2.3 we determine the bridging probability for the formation 
of a DNA bridge between two colloids, assuming the known thermodynamic 
parameters of DNA(hybridization free energy AG). Using this bridging probability 
as input, in section 2.4 wc calculate the effective binding free energy e^B for the 
formation of a dimer. Section 2.5 establishes the connection between the theory and 



^This restriction to binary systems is consistent with the current experimental approach. In a 
later chapter we will demonstrate that if each particle has a unique linker sequence, one might be 
able to programmably self-assemble nanoparticle clusters of desired geometry [65]. 
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Figure 2.1: Graphical depiction of various schemes for DNA bridging. A) A freely-jointed, rigid 
bridge constructed from complementary linker DNA. B) A flexible bridge can be constructed 
using complementary linker DNA. C) A rigid bridge constructed from short, flexible linker DNA 
and a long, rigid linker. 

the experimentally determined melting profile f{T), the fraction of unbound particles 
as a function of temperature. In particular, we demonstrate how knowledge of e^s 
can be used to determine this profile, including the effects of particle aggregation. 
In section 2.6 the theory is compared with two recent experiments detailing the 
reversible aggregation of colloids with DNA- mediated attraction [2], [1]. Section 
2.7 presents a detailed description of how the results can be applied to fit the 
experimental melting curves for a binary system of DNA-grafted colloids. The main 
results of the model are summarized in section 2.8. 

2.2 Description of the Problem 

We consider particles of type A and B which form reversible AB bonds as 
a result of DNA hybridization. The task at hand is to determine the relative 
abundance of the various colloidal structures that form as a function of temperature. 
From this information we can determine which factors affect the melting and 
aggregation properties in DNA-coUoidal assemblies. To do so we must determine 
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the binding free energy for all of the possible phases(monomer, dimer, infinite 
aggregate), and then apply the rules for thermodynamic equilibrium. As we will 
see, these binding free energies can all be simply related to cab, the binding free 
energy for the formation of a dimer. Our task is thus reduced to determining eab 
from the thermodynamic parameters of DNA and structural properties of the DNA 
linkers. In our statistical mechanical framework, cab is calculated from the model 
partition function, taking into account the appropriate ensemble averaging for the 
non-ergodic degrees of freedom. The result is related to the bridging probability 
for a pair of hnkers. By considering the specific properties of the DNA bridge that 
forms, the bridging probability can be related to the hybridization free energy AG 
of the DNA. In this way, we obtain a direct link between DNA thermodynamics 
and the global aggregation and melting properties in coUoidal-DNA systems. 

2.3 Bridging Probability 

To begin we relate the hybridization free energy AG for the DNA in solution to 

the bridging probability for a pair of linkers. This bridging probability is defined as 
the ratio ^"'^"^ , with Pbound the probability that the pair of linkers have hybridized 

-'jree 

to form a DNA bridge, and Pfree the probabihty that they are unbound. This ratio 
is directly related to the free energy difference of the bound and unbound states of 
the linkers AG (throughout this thesis we will use units with ks = 1): 



p 

-'■ bound 
free 



-AG 

T 



exp 



-AG 



(2.1) 



^ ^ /P(ri,r)P(r2,r)c^^r 

Here Co — IM is a reference concentration. P(r, r') is the probability distribution 
function for the linker chain which starts at r' and ends at r. The overlap density Cg// 
is a measure of the change in conformational entropy of the linker DNA as a result 
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of hybridization. It will depend on the properties of the linker DNA(ex: flexible vs. 
rigid), and the scheme for DNA bridging(ex: hybridization of complementary ends 
vs. hybridization mediated by an additional linker). Ce// is the concentration of free 
DNA which would have the same hybridization probability as the grafted linkers in 
our problem. As discussed in section 2.6, the DNA linker grafting density also plays 
an important role in determining the possible linker configurations and hence Ce//. 




Figure 2.2: The statistical weight of a bound state is calculated by determining the number of 
hybridized configurations for two complementary linker chains relative to the number of 
unhybridized configurations. 

Assuming that the size of the linkers is much smaller than the particle radius R, 
we first consider the problem in a planar geometry. Let the two linkers be attached 
to two parallel planar surfaces separated by a distance 2h. Referring to figure 2.2 we 
see that r' is the location where the linker DNA is grafted onto the particle surface, 
and r is the position of the free end. 

2.3.1 Hybridization Scheme A: Freely- Jointed Rigid Linkers 

In this section we consider hybridization by complementary, rigid linker DNA 
(scheme A in Figure 2.1). This scheme is particularly interesting since it is directly 
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related to several recent experiments [1], [2]. We assume that L < Ip and L <^ R, 
where Ip ~ 50nm is the persistence length of ds DNA and L is the ds linker DNA 
length. In this regime, the linker chains can be treated as rigid rods tethered on a 
planar surface. The interaction is assumed to be point-like, in which a small fraction 
A/L of the linker bases hybridize. 

We can calculate the overlap density by noting that the integral in Eq. (2.2) is 
proportional to the volume of intersection of two spherical shells (red and blue circles 
in Figure 2.3) : 

_ 2nrA _ Q {L - \r[ - r',\ /2) 

' {27rL^A)' ~ 2nL^\r[-v',\ ' ^'"^^ 



here A — A^/ sin/3 and r = y^L^ — |r^ — r'2f /4 (see notations in Figure 2.3). We 
have used the fact that cos/3/2 = jr'^^ — /2L. c^ff and the binding probability 
arc largest when the linkers are grafted right in front of each other, i.e. when 
\r[ — r'gl ~ 2h. By taking the hmit h L we arrive at the following result for the 
corresponding "bridging" free energy 

AG A ^ AG A + T log [AttL^Co] . (2.4) 

This free energy remains nearly constant for any pair of linkers, as long as they 
can be connected in principle, i.e. |r'j^ ~ < 2L. This limits the maximum lateral 
displacement of the linkers: r± < 2\Jl? — h?^ and therefore sets the elective 
cross-section of the interaction: 

a = Trri = 47r (L^ - /i^) (2.5) 
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Figure 2.3: Cross-sectional view of the hybridization of two complementary rigid linker DNA. The 
overlap density is calculated in a planar approximation to the particle surface. 

2.3.2 Hybridization Scheme B: Complementary Flexible 
Linkers 

We will now consider scenario B of figure 2.1, hybridization of complementary, 
flexible linker DNA. This situation can be realized in experiment by choosing linker 
DNA(ss or ds), provided the chain length L » Ip, the persistence length. We 
perform the calculation in a planar approximation to the particle surface, which 
implies the particle radius R » Rg, the radius of gyration of the linker chain. In 
scenario B, we must also take into account the entropic repulsion Grep of the linker 
DNA which arises as a result of confining the chains between planar surfaces. Since 
we are working with Gaussian chains, we can use the result of Dolan and Edwards 
[67]. Making the appropriate modification to equation 2.1 the binding probability 
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for a pair of linkers is the following^. 



exp 



-AG 



B 



T 



-eff 



exp 



■{AGB + 2Grep) 



(2.6) 



Defining x = ^ with the planar surfaces separated by a distance 2/i, the free energy 



of repulsion has the following asymptotic behavior. 



Grep = -7" log 



2x 




~ T 




k=l, k odd 
2 



log I 2x\l — ] + — , 



"24^2 
X < 1 



Grep = -7" log 



l-2^(-l)'=+iexp[-6xV] 



fe=i 



-riog[l-2exp(-6a;2)] x>l 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



The details of the calculation are given in Appendix 1. The final result gives the 



behavior of the binding probability exp 



-AGe 
T 



between complementary, flexible 



linkers as a function of x and the separation between grafting points Ar' = r'^^ — r'g. 



AGb ^ AGb + T 



AGb ^ AGb + T 



3 / Ar 
4 



t j + log (Rl^o) + log (f ) 

+31og(x) + ^ 



-2 log(x) + 3x^ 



9 / 3 
4 V T 



for X < 1 (2.11) 



for X > 1(2.12) 



Interpolating between the two regimes, we can see from the figure that the minimum 
free energy is at x < 1. 



"^Grep gives the free energy for a single linker with one end grafted on the planar surface, and the 
other end free. The binding probability contains a factor of 2Grep since each DNA bridge is made 
by joining 2 linkers. 
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Figure 2.4: ^^^(^rr — o) ^ function of x for flexible DNA bridges. The figure is normalized by 



choosing + log(i?^Co) = 0. 



2.3.3 Hybridization Scheme C: Short Flexible Markers with 
a Long Rigid Linker 

We now turn our attention to scenario C of figure 2.1, hybridization of short, 
fiexible marker DNA with radius of gyration Rg to a long, rigid hnker DNA of length 
L. We will consider the case ^ = ^ 1. For this reason we can neglect the 
entropic repulsion Grep of the short linkers, since they only feel the presence of the 
surface to which they are attached. However, in this scenario we must take into 
account the loss of entropy of the long, rigid DNA linker. After hybridization this 
linker strand does not have the full 47r steradians of rotational freedom it does when 
free in solution. The appropriate modification to the binding probability is: 



exp 



-AG, 
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'eff 



exp 
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(2.13) 



(2.14) 
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Once again, the reader interested in the details of the calculation is directed to 
Appendix B. For completeness we quote the result here. 



2.4 Effective Binding Free Energy 

We now proceed with the calculation of the effective free energy eAB , which 
is associated with the formation of a dimer from a pair of free particles, A and 
B. Since the DNA coverage on the particle surface is not uniform, this free 
energy, and the corresponding partition function Z, would in principle depend on 
the orientations of the particles with respect to the line connecting their centers. 
The equilibrium binding free energy would correspond to the canonical ensemble of 
all possible orientations, i.e. cab — —T\ogATr{Z). However, this equilibrium can 
only be achieved after a very long time, when the particle pair samples all possible 
binding configurations, or at least their representative subset. The real situation is 
different. After the first DNA- mediated bridge is created the particle pair can still 
explore the configurational space by rotating about this contact point. However, 
after the formation of two or more DNA bridges (at certain relative orientation of 
the particles) , further exploration requires multiple breaking and reconnecting of the 
DNA links, which is a very slow process. We conclude that the system is ergodic 
with respect to the various conformations of the linker DNA for fixed orientations of 
the particles, but the orientations themselves are non-ergodic variables. The only 
exceptions are the single-bridge states: the system quickly relaxes to a more favorable 
orientational state (unless the DNA coverage is extremely low, and finding a second 
contact is very hard). If N denotes the number of DNA bridges constituting the 
AB bond, the appropriate expression for €ab in this partially ergodic regime is the 



AG'c(Ar' 



0,e = e*) ~ AGc + Tlog 
AGc + 4.24r + Tlog [L^RgCo] 




(2.15) 



(2.16) 
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so-called component averaged free energy [68], [69]: 



CAB ^-T {log Z) 



N>2 



(2.17) 



Each DNA bridge between particles can be either open or closed. 



Zbridge = 1 + exp 



T 



(2.18) 



Here r ■ is the 2D position where the bridge is grafted onto surface i. We now consider 
a generic case when the interaction free energy AG depends on the separation 
between planar surfaces 2h', and the separation of grafting points r'^ — r'g, without 
assumption of a particular bridging scheme. If the probability for bridge formation 
is small, two DNA linkers on the same surface will not compete for complementary 
linkers. In this regime the free energy can be calculated by summing over the 
contribution from each bridge that forms between dimers. 

AG{h',r',-r') 



1 + cxp 



(2.19) 



We convert the summation to integration by introducing the linker areal grafting 
density a. 

AG(/i',r;-ry 




1 + exp 



T 



d%d% (2.20) 



Changing variables to Ar = r'^^ — r2 and p = (r'^ + /2, we can reintroduce the 
notion of a bridging cross-section a{h'), this time in a model- independent manner: 



a(h')\og 



1 -I- exp 



d^Arlog 



1 -I- exp 



AG{h', Ar) 
f 



(2.21) 



Here AGo{h') = AG{h', Ar = 0) is the minimum free energy with respect to the 
separation between grafting points Ar. We can now write the free energy: 

AGo{h' 



^^-^J <7i(p)^T2(p)a(/i')log 



1 -I- exp 



(2.22) 
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We now convert from the planar geometry to the spherical particle geometry using 
the Derjaguin approximation [70]. 

d^p = pdpd(j) (2.23) 
h'^h+^ (2.24) 

Let AG* be the minimal value of the bridging free energy. Then the result for F 
can be rewritten as: 



F^-TN log 



1 + exp 



AG, 
T 



(2.25) 



Here N has a physical meaning as the number of potential bridges for given relative 
positions and orientations of the particles: 



N = J ai{'p)a2Cp)a{h' 



log 



l + exp -AGo(/i')/T 



(2.26) 



log [l + exp (^-AG^/r) 

One can calculate the average value of in terms of the average grafting density, 
(7 = ((ji) (da) : 



(N) = 2nRa' 



■ ja{h') 



log l + exp -AG„(/i')/T 



dh' 



log 1 + exp (^-AG*/r) 

In a generic case of randomly grafted linkers, (A^) completely defines the overall 



(2.27) 



distribution function of A^, which must have a Poisson form: P{N) 



The average number of bridges {N) between two particles depends on both the DNA 
linker grafting density a and the bridging probability determined from AG. 

The free energy for the formation of a dimer eab — {F)2+ ~ T log fl. The second 
term is the entropic contribution to the free energy, which comes from integration 
over the oricntational and translational degrees of freedom of the second particle. 
Because the system is not ergodic in these degrees of freedom, the accessible phase 
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space Q will be reduced by a factor of P2+. P2+ is the probability that there are at 
least two DNA bridges between the particles. In terms of the average number of 
bridges {N) between particles, we have the following relations: 



P2+ = 1 - (1 + (iV))e-W 
{N) (1 - e-W) 



eAB^-T{{N)^. log 



1 + exp 



P2^ 

T 



(2.28) 
(2.29) 



+ \og[P2^AT:5{2Rfco] ) (2.30) 



Here 5 is the localization length of the AB bond, which comes from integrating the 
partition function over the radial distance between particles. 

2.4.1 Scheme A 

We now can calculate {N) for the case of freely-jointed rigid bridging considered 
earlier (i.e. for scheme A). In a previous section we provided a direct calculation 
of the interaction free energy, AGo{h) ~ const = AG^ (eq.2.4), and bridging 
cross-section, a{h') — 47r(L^ — h''^). Applying eq. 2.27 we arrive immediately at the 
following result. 

L 

(N) = STr'a'R J {L^ - h")dh' = ^^vrV^i^L^ ^^^^^ 



2.4.2 Scheme B 

We note that in this case, since the binding probability for a given pair of linkers 
is Gaussian in the separation between grafting points Ar = r'^^ — we can perform 
an analytic calculation of the effective cross section. In what follows x' — h'/Rg. 
Recall the definition of a{h'): 



a{h')log 
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(2.32) 



27 



exp 



T 



A{x') exp 



3 /Ar 

4 



-(AGs + '^Grep) 
f 



(2.33) 



(2.34) 



The explicit form of the dimensionless concentration c{x') is given in Appendix A 
(see Eq. A. 11). Changing to polar coordinates (r, 9) we have: 

2tt oo 



a(a;') log [1 + A(x')] = J d9 J rlog 







1 + A{x') exp 



3 / r 
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dr (2.35) 



Define a new variable u{r) = A{x') exp 



3 r 



4 \ R. 



log [1 + u] 4:71 



a(x')log[l + A(x')] = -yi?^ J du^^^^^;^^-^RlLi,[-A{x')] (2.36) 



A(a;') 



The calculation yields the following result, with Li2 [z] = |^ the Dilogarithm. 



fc=i 



a{x') 



An Lt2 [-A{x')] 



(2.37) 



3 log [1 + A{x')] ^ 

Since Li2{z) < for z < 0, a{x') is positive as required. Prom this effective cross 
section we can compute the average free energy {F) as a result of DNA bridging 
between particles, with a the average areal grafting density of DNA hnkers. Here 
p — r[ + r'2, with r[ the location where hnker i is grafted on the planar surface. 



(F) = -Ta' J a{x') log [l + A(x')] d'^p 



(2.38) 



Converting from the planar geometry to the spherical nanoparticle geometry using 
the Derjaguin approximation we have: 

^2 



X = X + 



IRRn 



(2.39) 
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CXD 

.2 



(F) = T^a^RRl I Li2 [-A(a;')] dx (2.40) 



This integration can be performed numerically. As discussed previously, the 
free energy cab for the formation of a dim.er{AB pair) also contains an entropic 
contribution from integration over the orientational and translational degrees of 
freedom of the second particle. 

2.4.3 Scheme C 

We can also determine the free energy in this scenario using the approximation 
method developed. 

AGc( Ar' = 0, e = e* 



F^-T (N) log 



1 + exp 



T 



(2.41) 



We provide a simple geometrical argument to determine the average number of DNA 
bridges (N) between particles. We assume that the rigid hnkers are aligned with a 
small component parallel to the surface. 

a{h) ~ V (2.42) 
y = (L + A)tane^ax~i?g(l + e) (2.43) 



Then applying equation 2.27 with h — ^(1 + e) we have: 

(iV) = n^a^RR^L J {1 + tfde ~ 2V27rV'i?i?^ (2.44) 



2.5 Aggregation and Melting Behavior 

At this stage we have calculated the binding free energy eab for an AB pair, 
starting with the thermodynamic parameters of DNA (hybridization free energy 
AG). In this section we establish the connection between that result and the 
experimentally observable morphological behavior of a large system. One of the 
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ways to characterize the system is to study its melting profile f{T), which is the 
fraction of unbound particles as a function of temperature. To determine the profile 
we calculate the chemical potential for each phase(monomer, dimer, etc.) and apply 
the thermodynamic rules for phase equilibrium. We will demonstrate how the single 
binding free energy cab can be used to determine the contribution of each phase to 
the melting profile, including the effects of aggregation. 

2.5.1 Dimer Formation 

To begin we discuss the formation of dimers via the reaction A + B AB. 
We can express the chemical potential of the i*'* species //^ in terms of the particle 
concentrations Ci — ^. 

fiA^T\og{cA) (2.45) 
l^s log (cb) (2.46) 
II AB = Tlog (cab) + eAB (2.47) 

Here sab is the binding free energy for the formation of a dimer. In terms of the 
potential V{r) between A and B type particles we have: 

V{r 



eAB = -T log 



4:Ti{2Rfco J drexp 



(2.48) 



T 

In this section we are not particularly concerned with the specific form of the 
DNA-induced potential V{r), having already determined cab in the previous section. 
We simply note that the prefactor 47r(2i?)^ arises since the interaction is assumed to 
be isotropic, with R the particle radius. Equilibrating the chemical potential of the 
various particle species, we obtain the condition for chemical equilibrium. 

I^A + I^B^ I^AB (2-49) 
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The result is a relationship between the concentration of dimers and monomers. 



caCb 
Cab = exp 



(2.50) 



The overall concentration of particles in monomers and dimers must not differ from 
the initial concentration. 



c\ = ca + Cab 



cb + Cab 



(2.51) 
(2.52) 



If the system is prepared at equal concentration, c\ — — \ctot-i subtracting the 
two equations we see that ca = cb = c. Written in terms of the fraction of unbound 



particles / 



we have a quadratic equation for the unbound fraction. 

2 Ctot 



1 = / + exp 



—^AB 

T 



(2.53) 



To simplify we have defined an effective free energy cab for the formation of a 
dimer. 



€aB = eAB-T log 



Ctot 

2^ 



(2.54) 



The solution for the fraction of unbound particles as a function of temperature is 
simply: 



1 + -^1 +4exp [^] 
2 exp 



(2.55) 



Previous studies [1] only included the dimer contribution to the melting properties 
of DNA colloidal assemblies. With the basic formalism at hand, we can now extend 
the preceding analysis to include the contribution of trimers and tetramers. 
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2.5.2 Trimers and Tetramers 

Now consider the formation of a trimer via 2A + B ABA. The chemical 
potential is shghtly different in this case. 

f^ABA = Tlog (caba) + €.ABA (2.56) 

Taking into account that there are now two AB bonds in the structure, one might 
conclude that caba = '^^ab- This is not quite correct, since there is a reduction in 
solid angle available to the third particle. To form a trimer, an AB bond forms first, 
which contributes eAs to eASA- Some simple geometry shows that the remaining A 
particle only has Stt steradians of possible bonding sites to particle B. Making this 
change in the prefactor of eq. 2.48, one can see that the second bond contributes 
eAB - T\og (I) to eABA- 

eABA = '^eAB -T log (2.57) 

The equation for chemical equilibrium can once again be expressed in terms of the 
particle concentrations. 

+ /^s = ^^ABA (2-58) 



3 c^cb 
CAB A = 7— ^-exp 



-2eAi 
T 



(2.59) 



4 4 

To include the trimer contribution, we note that there are two possible varieties, 

with eABA = ^BAB- 



c\^ca + Cab + '^caba + cbab (2.60) 
Cb = cb + Cab + caba + '2cbab (2.61) 
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Following the same line of reasoning as before, the resulting equation for the unbound 
fraction / is: 



1 = / + exp 



f + 7exp 



-2? 



AB 



T 



f 



(2.62) 



For tetramers we will follow the same general reasoning, however in this case 
there are two different structure types. The reaction 2A-\-2B ^ ABAB results in 
the formation of string like structures. 



I^ABAB — T\og {cabAb) + ^ABAB 



(2.63) 



As in the trimer case, the last particle has 37r steradians of possible bonding sites, 
and contributes cab — T\og (|) to eASAB- 

2" 



^ABAB = 3eAB - T\og 



2fJ,^ + 2/^5 — l^ABAB 
CABAB = 



T 



(2.64) 

(2.65) 
(2.66) 



If an A type particle approaches a trimer of variety ABA, a branched structure can 
result. The reaction 3 A + B ^ AAAB results in the formation of these branched 
structures. 

I^AAAB = Tlog (caaAb) + (^AAAB (2.67) 

For the branched case, the last particle has approximately 27r steradians of possible 
bonding sites, and contributes cab — Tlog (|) to €aaab- 

(2.68) 
(2.69) 

caaab = I - I exp — (2.70) 



^AAAB = 3eAB ~ 
SHa + l^B — t^AAAB 



3\ CaCb 



-3eAi 



To include all of the tetramer contributions, note that there are two branched 
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varieties, with caaab — ^bbba- Finally we impose the constraint that the initial 
particle concentrations do not differ from the concentration of all the n-mers, for 
n=l,2,3,4. 

c\ — ca-\- Cab + 2caba + cbab + '^cabab + Sc^aas + cbbba (2-71) 
Cfi = cs + Cab + Cab A + '^cbab + '^cabab + caaab + "^cbbba (2.72) 

The final result is an equation for the unbound fraction / expressed entirely in terms 
of the effective free energy Iab of a dimer. 



1 = / + exp 



-^AB 



/' + ?exp 



-2?. 



AB 



f + 



21 



exp 



-3?, 



AB 



(2.73) 



For high temperatures, the melting profile is governed by the solution to this 
polynomial equation for /. For temperatures below the melting point we expect 
to find particles in large extended clusters. We now proceed to calculate the 
equilibrium condition between monomers in solution and the aggregate. 

2.5.3 Reversible Sol-Gel Transition 

To understand the basic structure of the aggregate, we simply note that there 
are many DNA attached to each particle. This gives rise to branching, as in the 
discussion of possible tetramer structures. Since the DNA which mediate the 
interaction are grafted onto the particle surface, once two particles are bound, the 
relative orientation of the pair is essentially fixed. The resulting aggregate is a 
tree-like structure, and the transition to an infinite aggregate at low temperatures is 
analogous to the sol-gel transition in branched polymers [21]. 

Particles in the aggregate are pinned down by their nearest neighbor bonds, so 
we do not consider their translational entropy. As a result the chemical potential is 
simply = €00- Equilibrating the chemical potential of the monomer in solution 
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Figure 2.5: The actual unbound fraction / is the concatenation of the aggregate profile for 
T < T* and the n-mer profile for T > T* . The fraction of particles in dimers, trimers, and 
tetramers is also plotted. 

and in the aggregate we have: 

Tlog(c)=eoo (2.74) 
Coo = eAB - Tlog (7^) (2.75) 

Here 7^^ ~ 1 is the configurational entropy of the branched aggregate, per particle. 

The concentration of particles in the aggregate Coo is the the total concentration 
minus the n-mer concentration. Here Ci = ca + Cb is the total monomer 
concentration, C2 = cab is the total dimer concentration, etc. 

Coo ~ Ctot - Ci - C2 - C3 - C4 (2.76) 

Expressed in terms oi'eAs and the fraction of solid angle available to particles in the 
aggregate = ^ we have: 

(2.77) 



, 1 

Joo = — exp 

O'oo 
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The transition from dimers, trimers, etc. to the aggregation behavior is the 
temperature T* at which f^{T*) is a solution to eq. 2.73. In words, T* is the 
temperature at which the aggregate has a non-zero volume fraction. The fraction 
of unbound particles for these colloidal assembhes will be governed by eq. 2.77 for 
T < T* and eq. 2.73 for T > T*. As claimed, we can simply relate the unbound 
fraction to e^^ for both n-mers and the aggregate. 

2.6 Comparison to the Experiments 

Let's consider the experimental scheme of Chaikin et al [1]. In the experiment, 
R = .5fim polystyrene beads were grafted with ds DNA linkers of length L ~ 20nm. 
The 11 end bases of the A and B type particles were single stranded and 
complementary. We have already determined the bridging probability in this 
scenario(see scheme A). In the experiment [1] a polymer brush is also grafted onto 
the particle surface, which will have the effect of preferentially orienting the rods 
normal to the surface(See Figure 2.3). This confinement of the linker DNA can be 
incorporated quite easily into our results for AG and (N). To modify Eq. 2.4, when 
integrating over linker conformations we simply confine each rigid rod to a cone of 
opening angle 2a. The upper bound for the polar integration is now a as opposed 

to TT. 



The alignment effect should also be taken into account when calculating (N). If 
the particles are separated by less than 2L cos a the end sequences will be unable to 



AG A ^ AG A + Tlog [4nL^Co{l - cosaf] 



(2.78) 



hybridize. Following the same steps as before, the lower bound for the h' integration 



is now Lcosa as opposed to 0. 




3)] 



(2.80) 



(2.79) 
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In the absence of the brush, and at sufficiently low linker grafting density a, the 
alignment effect could be removed by setting a = |, in which case we recover 
our previous results. Since the polymer brush is stiff, it also imposes a minimum 
separation of 2h between particles, where h is the height of the brush. As a result, 
in the expression for eAB we can approximate the radial flexibility of the AB bond 
as 5 ^ L — /i. 




TemperaturefC) 

Figure 2.6: Comparison of the melting curves /(T) determined by our model to the experimental 
data of Chaikin et al(See Fig. 2 in [1]). The four data sets are for the four different polymer 
brushes used. For the model fits we find that (A^)2^=2.01 for crosses, 2.07 for solid triangles, 2.f3 
for empty triangles, and 2.35 for squares. 

We have now related the free energy eAB to the known thermodynamic 
parameters of DNA(AG' = AH - TAS, AH = -77.2^ and AS = -227.8^), 
and the properties of linker DNA chains attached to the particles (grafting density 
cr ~ 3 X 10^ ^^^ and linker length L ^ 20nm). The height of the polymer brush 
is = 13 ± 5nm [1]. In fitting the experimental data we have taken the average 
value (h) = 13nm. Changing h within these bounds does not have a major effect 
on the melting curves. As a result there is one free parameter in the model, the 
confinement angle a. This angle determines (A^) and AG, which in turn determine 
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e^AB, and finally the melting profile /. 

With some minor modifications we can also analyze the "tail to tail" hybridization 
mode in a recent experiment of Mirkin et al [2]. In this experiment, R = 6.5nm 
gold nanoparticles were chemically functionalized with ss DNA linkers. The last 15 
bases on the markers for particles of type A and B were chosen to be complementary 
to a 30 base ss DNA linker. Since the strands are not ligated after hybridization, 
the experimental pictures are similar. 




TemperaturefC) 

Figure 2.7: The effect of the linker DNA grafting density a on the melting profile /(T). The 
results of the model are compared with experimental data in [2] . The three data sets represent 
grafting densities of 100% (squares), 50%(circles), and 33%(triangles) for which (A^)2^=2.32, 2.16, 
and 2.05 respectively. 

The unhybridized portion of the ss DNA linker simply serves as a spacer, and the 
hybridized portions become ds DNA, which we can again treat as rigid rods. This 
experiment is done without the addition of a polymer brush, but the grafting density 
is two orders of magnitude larger than the experiment of Chaikin et al. As a result, 
there is still an entropic repulsion [55] associated with compressing the particles 
below separation 2h. Here h could loosely be interpreted as the radius of gyration 
of the unhybridized portion of the linker. Despite the fact that L ~ i?, our planar 
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calculation of AG provides a good fit to the experimental data. The other major 
difference is that now the attraction between particles is mediated by an additional 
DNA linker. 



AG^ AGi + AG2-T\og 



(2.81) 



The term AGi(AG'2) is the contribution to the free energy from the hybridization 
of the linker on an A{B) type particle to the complementary portion of the 30 
base ss linker. The hybridization free energies AGi and AG2 were calculated 
with the DINAMelt web server [71]. The last term is the contribution to the free 
energy from the translational entropy of the additional linker DNA, with cunk the 
additional linker concentration. This highlights some incorrect assumptions of the 
thermodynamic melting model [2], where the two hybridization free energies were 
not calculated separately, and the translational entropy of the additional linker DNA 
was ignored. By introducing dilutent strands to the system, one can probe the effect 
of the hnker grafting density a on the melting properties of the assembly(See Figure 
2B in [2]). The agreement between the experimental data and our theory is good, 
except at small / values. This is not surprising, since comparing the two requires 
relating the measurement of optical extinction to the unbound fraction /. This is a 
nontrivial matter when dealing with aggregation, which corresponds to the small / 
regime. 

2.7 Fitting Algorithm 

In this section we present a step by step method for fitting the melting curves 
obtained experimentally for a binary system of DNA-grafted colloids. 
Step 1: Determine AG 

The first step is to determine the hybridization free energy AG for the 
DNA strands free in solution. In many cases the value has been determined 
experimentally. Alternatively, there are a number of web based applications which 
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calculate hybridization free energies. For example, the DINAMelt server which 
can be located at http://www.bioinfo.rpi.edu/applications/hybrid/hybrid2.php 
and NUpack which can be located at http://piercelab.caltech.edu/nupack. Note 
that in the case where the hybridization is mediated by an additional linker, the 
translational entropy of that linker must be taken into account (see Eq. 2.81). 
Step 2: AC ^ AG 

Since the DNA linkers in our problem are grafted onto the particle surface, 
we need to determine how the grafting effects the bridging probability (see Eqs. 
2.1 and 2.2). This entails calculating the overlap density Cg// which is a measure 
of the change in conformational entropy of the DNA strands upon hybridization. 
Determining the appropriate calculation will depend on the hybridization scheme 
(see Fig. 2.1). In this chapter calculations have been performed for three different 
schemes (see Eqs. 2.4, 2.11, 2.15), although the effects of linker confinement have 
only been taken into account in scheme A (see Eq. 2.78). 

Step 3: Calculate (N) 

The next step in the procedure is to determine the average number of bridges 
(A^) that form between an AB pair. The general starting point is Eq. 2.27. In this 
chapter calculations have been performed for three different hybridization schemes 
(see Eqs. 2.31, 2.40, 2.44). The effects of linker confinement have been taken into 
account in scheme A (see Eq. 2.79). Note that in our approximation scheme the 
general relation between (N) and the free energy F is given by Eq. 2.25. 

Step 4: Determine €ab 

The next step in the procedure is to relate (N) to the binding energy for the 
formation of a dimer pair eab (see Eq. 2.30). The quantity of interest for the fitting 
e^B is simply related to €ab by Eq. 2.54. 

Step 5: Determine the melting profile /(T) 

We are now in a position to relate the calculation to the experimentally measured 
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quantity f(T), which is the fraction of monomers as a function of temperature. For 
high temperatures / is determined by the solution of the polynomial Eq. 2.73. As 
the temperature is lowered at T = T* we reach the point where / determined from 
the n-mer profile (Eq. 2.73) is equal to / determined from the aggregate profile (Eq. 
2.77). For T < T* the melting profile is determined by Eq. 2.77. 

2.8 Summary 

We have developed a statistical mechanical description of aggregation and 
melting in DNA-mediated colloidal systems. First we obtained a general result for 
two-particle binding energy in terms of DNA hybridization free energy AG, and 
two model-dependent parameters: the average number of available bridges {N) and 
the overlap density for the DNA Cg//. We have also shown how these parameters 
can be calculated for a particular bridging scheme. In our discussion we have 
explicitly taken into account the partial ergodicity of the problem related to slow 
binding-unbinding dynamics. 

In the second part it was demonstrated that the fractions of dimers, trimers and 
other clusters, including the infinite aggregate, are universal functions of a parameter 
'^ab/T — €ab/T — log [ctot/2co]. The theory has been calculated for three separate 
hybridization schemes. The obtained melting curves are in excellent agreement 
with two types of experiments, done with particles of nanometer and micron sizes. 
Furthermore, our analysis of the experimental data give an additional insight into 
microscopic physics of DNA bridging in these systems: it was shown that the 
experiments cannot be explained without the introduction of angular localization of 
hnker dsDNA. The corresponding locahzation angle a is the only fitting parameter 
of the model, which allows one to fit both the position and width of the observed 
melting curves. 

There are several manifestations of the greater predictive power of our statistical 
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mechanics approach, compared to the earher more phenomenological models. First, 
once a is determined for a particular system, our theory allows one to calculate 
the melting behavior for an alternative choice of DNA linker sequences. Second, 
if the resulting clusters are separated, for example in a density gradient tube, the 
relative abundance of dimers, trimers, and tetramers can be compared to the values 
determined from the theory. 

Finally, the theory predicts aging of the colloidal structures, one experimental 
signature for which is hysteresis of the melting curves. Such an experiment proceeds 
by preparing a system above the melting temperature, and measuring the unbound 
fraction of colloids as the temperature is lowered. The system is allowed to remain in 
this cooled state for a very long time, perhaps months, during which multiple DNA 
bridges break and reform. During this time the colloids relax into a more favorable 
orientation state, including states which are not accessible by simply rotating about 
the contact point formed by the first DNA bridge between particles. This favorable 
orientation state is characterized by an average number of DNA bridges {N) greater 
than what we calculate in the partially ergodic regime. If the unbound fraction 
is then measured as the temperature is increased, the melting curve will shift to a 
higher temperature, consistent with a larger value of (A^") . 
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CHAPTER III 



DYNAMICS OF "KEY-LOCK" INTERACTING 

PARTICLES 

3.1 Introduction 

In this chapter [72] , [73] we present a theoretical study of desorption and diffusion 
of particles which interact through key-lock binding of attached biomolecules. It is 
becoming common practice to functionalize colloidal particles with single-stranded 
DNA (ssDNA) to achieve specific, controllable interactions [26], [1], [43], [15], [25], 
[49]. Beyond the conceptual interest as a model system to study glassiness [65] and 
crystallization, there are a number of practical applications. Colloidal self-assembly 
may provide a fabrication technique for photonic band gap materials [74], [75]. 
One of the major experimental goals in this line of research is the self-assembly 
of colloidal crystals using DNA mediated interactions. The difficulty stems in 
part from the slow relaxation dynamics in these systems. The main goal of this 
chapter is to understand how the collective character of key-lock binding influences 
the particle dynamics. In doing so we gain valuable insight into the relaxation 
dynamics, and propose a modified experimental setup whose fast relaxation should 
facilitate colloidal crystallization. 

Similar systems have also attracted substantial attention in other areas of 
nanoscience. In particular, by functionalizing nanoparticles with antibodies to a 
particular protein, the nanoparticles have potential applications as smart, cell-specific 
drug delivery vehicles [76], [77]. These nanodevices take advantage of the fact 
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that certain cancerous cells overexpress cell membrane proteins, for example the 
folate receptor. An improved understanding of desorption and diffusion on the cell 
membrane surface may have implications for optimizing the design of these drug 
dehvery vehicles. This is the subject of chapter 7. 

In what follows we present our results on the dynamics of particles which interact 
through reversible key-lock binding. The plan for the chapter is the following. In 
section 3.2 we introduce the key- lock model and explain the origin of the two model 
parameters A and m. The parameter A determines the binding energy for the 
formation of a key-lock pair. The parameter m is the mean of the distribution 
for the number of key-lock bridges. Depending on m, which is related to the 
coverage of the functional groups (e.g. ssDNA), there are two distinct regimes. 
At low coverage there is an exponential distribution of departure times, but no 
true lateral diffusion. As the coverage increases, we enter a regime where the 
particle dynamics is a result of the interplay between desorption and diffusion. An 
estimate is provided for the value of m which determines the crossover from the 
localized to diffusive regime in section 3.3. In section 3.4 the localized regime is 
discussed in detail. In this regime the particle is attached to a finite cluster and 
remains localized near its original location until departing. We derive the partition 
function for the finite clusters, and calculate the departure time distribution. In 
section 3.5 we determine the departure time distribution in the diffusive regime. 
We present an effective Arrhenius approximation for the hopping process and a 
Fourier transform method which greatly simplifies the calculation. In section 3.6 
we discuss the random walk statistics for the particles' in-plane diffusion. A set 
of parametric equations is derived to relate the average diffusion time to the mean 
squared displacement. The lateral motion is analogous to dispersive transport in 
disordered semiconductors, ranging from standard diffusion with a renormalized 
diffusion coefficient to anomalous, subdiffusive behavior. In section 3.7 we connect 
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Figure 3.1: Graphical depiction of particles interacting with a flat, 2D substrate by multiple 
key-lock binding. 

our results to recent experiments with DNA-grafted colloids. We then discuss the 
implications of the work for designing an experiment which facilitates faster colloidal 
crystallization. In section 3.8 we conclude by summarizing our main results. 

3.2 Model Description 

We now present the model, where a single particle interacts with a flat two- 
dimensional surface by multiple key lock binding (see Fig. 3.1). At each location 
on the surface there are m key-lock bridges which may be open or closed, with a 
binding energy of e for each key-lock pair. Here we have neglected the variation in 
e. In the case of the DNA-coUoidal system mentioned in the introduction, the model 
parameter e is related to the hybridization free energy of the DNA. The resulting 
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m-bridge free energy plays the role of an effective local potential for the particle [30] : 



Generically, m is a Poisson distributed random number — m'"exp(— m)/m! 
where m denotes the mean of the distribution. The model parameter m is a 
collective property of the particle-surface system. For example, consider the case 
of dendrimers functionalized with folic acid, which can be utilized for targeted, cell 
specific chemotherapy. The folic acid on the dendrimer branch ends form key-lock 
bridges with folate receptors in the cell-membrane. In this case rn will depend on 
the distribution of keys (folic acids) on the dendrimer, and the surface coverage of 
locks (folate receptors) in the cell membrane. 

At each location, the particle is attached to the surface by m bridges. To detach 
from the surface the particle must break all its connections, in which case it departs 
and diffuses away into solution. Alternatively the particle can hop a distance a to a 
new location characterized by a new value of the bridge number m. By introducing 
the correlation length a, we have coarse-grained the particle motion by the distance 
after which the new value of the bridge number becomes statistically independent of 
the value at the previous location. In the localized regime the particle remains close 
to its original location until departing. In the diffusive regime the particle is able to 
fully explore the surface through a random walk by multiple breaking and reforming 
of bridges. 



U{m) = -TmA 



(3.1) 



A = log(l + exp[e/r]) 



(3.2) 
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Figure 3.2: The crossover from the locahzed to the diffusive regime below the percolation 
threshold. Estimates based on the characteristic cluster size (Eq. 3.3, red line) and confinement of 
the random walk (Eq. 3.4, blue line) give similar crossovers. For large A the crossover condition is 



3.3 Crossover from Localized to Diffusive Behav- 



Naively one might expect the crossover between the two regimes to occur at 
the percolation threshold, where one first encounters an infinitely connected cluster 
of sites with m > 0. However, the crossover from the localized to diffusive regime 
occurs at smaller m than predicted by percolation theory. If pc = 1/2 denotes 
the critical probability for site percolation on the triangular lattice, the percolation 
transition occurs at m = log (2). There are two alternative estimates for the 
crossover from the localized to the diffusive regime. The first is to compare the 
average number of steps n = exp (Am) the particle takes before departing (see 
section 3.5) to the characteristic cluster size Sc = l/log(l/Ap) below the percolation 
threshold. Here A = 5.19 is a numerical constant for the triangular lattice [78], and 
in the percolation language p = 1 — exp(— m) is the occupancy probability. The 
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crossover condition n — Sc can be expressed as a function of m. 

^--^^og[-log{X{l-e-^)}] (3.3) 

Alternatively, in the localized regime the particles' random walk is confined by the 
characteristic cluster size. Below percolation the radius of gyration of the cluster is 
Rs ~ with p — 0.641 in two dimensions. Comparing the radius of gyration of the 
cluster to the radius of gyration for the particles' random walk, the crossover occurs 
at n = (sc)^''- 

A = -|log[-log{A(l-e--)}] (3.4) 

fit 

Since 2p differs from 1 by less than 30%, both conditions give similar crossovers (see 
Fig. 3.2). The saturation at m = — log (l — j) occurs for very large A, as a result 
for binding energies of a few T per bridge the crossover occurs at m ~ 0.1. 

3.4 Localized Regime 

In the percolation language, when the occupancy probability p = 1 — exp(— m) is 
small, particles are localized on finite clusters. In this localized regime particles are 
able to fully explore the cluster to which they are attached before departing. This 
thermalization of particles with finite clusters permits an equilibrium calculation 
of the cluster free energy F — —Tlog{Z). The departure rate is given by the 
Arrhenius relation K — ■^exp{F/T). Here Tq is a characteristic timescale for 
bridge formation. The probability that the particle departs between t and t + dt is 
determined from the departure time distribution ^{t)dt ~ K exp[—Kt]dt. 

To begin we calculate the partition function for the finite clusters. The cluster 
is defined as s connected sites on the lattice, all of which are characterized by 
< m < m* bridges. For Poisson distributed bridge numbers the partition function 
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for the finite cluster is: 

s m* — l 



Z(m*,s) = Pm,exp{Ami) = — — — — -(exp(me^)Q(me^, m*) - 1) (3.5) 



s 

Prn, exp(Ami) = 7 

eifnl 

1=1 mi=l 

Because by definition the cluster does not contain sites with m = bridges we 
have renormalized the probability distribution Pm = Pm/{^ — exp(— m)) so that 
E^=i Pm = 1- Here Q{x, m*) = r{x, m*)/r{m*) = exp(-x) ^^=0"' ^V^' is the 
regularized upper incomplete F function. In the language of the statistics of extreme 
events, m* — 1 is the maximum "expected" value of m in a sample of s independent 
reahzations [79] . The point is that on finite clusters we should not expect to achieve 
arbitrarily large values of the bridge number. Hence when averaging the partition 
function to obtain the cluster free energy one should only average over sites with 
m < m*. The distribution function for m* is obtained by noting that the probability 
that all s values of m are less than m* is {^Zl'i' Pm)' = ( ''""T^^p'jg-f By 
differentiating this quantity with respect to m* we obtain the distribution function 
for the maximum expected value of m. 

f( *N /exp(m)g(m,m*) - iV"^ ~ 

fs{m = s J— Pm- 3.6) 

\ exp(m) — 1 / 

The cluster size distribution below the percolation threshold is exponential [80] with 
characteristic cluster size Sc — 1/ log(l/Ap). 

P.(p) = ^^exp(^-^) (3.7) 

The summation over s can be performed analytically, which allows the result to be 

expressed as a single summation over m* . 

00 00 

m*=2 s=l 

A(l - A(l - exp(-^))) ^ p^.(,^(^,.)Q(^,A_„.) _ 1) 1 + 

exp(m) - 1 ^2 (1 - y{m*)f 
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Figure 3.3: Departure time distribution function versus time in the localized regime with m — 0.1. 
The results of our calculation (solid lines) are compared to single exponential relaxation with 
departure rate K — "^p^"'^) (dotted hues). 

y{m*) = X{Q(jn,m*) — exp{—m)) (3.9) 

For fixed fn, as in the plot (see Fig. 3.3), changing A is directly related to a 
change in the average binding free energy. Increasing A leads to a reduction in the 
rate of particle departure. 

3.5 Diffusive Regime 

The departure time distribution changes significantly in the diffusive regime. In 
this regime the particle can explore the surface to find a more favorable connection 
site, which leads to a longer lifetime for the bound state. This phenomenon is 
qualitatively similar to aging in glassy systems. In these systems one finds that 
the response to an external field is time dependent [81]. In the magnetic analogy 
this leads to a time dependence of the magnetization. Below the glass temperature, 
the longer one waits before applying the external magnetic field, the more time the 



50 



10 



15 



10 



10 



Tr 



m 



TO 



10^ 



10 



. Triangular Lattice 






m = 1 / , 








m = 2 ^ 


• 


X * 




m =4 

^^^^^^^ 
^^^^^^^ 


X • 

M * 


^^^^^^^ 

^ • 


• 

m =8 









10 



15 



Figure 3.4: Comparison of the ensemble averaged dwell time in a lattice model (Eq. 3.10) to the 
effective Arrhenius approximation (Eq. 3.11). In the plot the product Am = 4 is held constant. 



system has to settle into deep energy wells, and the smaller the response. In our 
case, the diffusive exploration of the particle allows it to find a deeper energy well, 
which leads to an increase in the bound state lifetime. 

As a result, the departure time distribution must now reflect not only desorption, 
but also hopping to adjacent sites. The hopping rate between neighboring sites % 
and j is given by an Arrhenius law = :^exp[— A(mj — m^QivHi — mj)], with 
Q[x) the Heaviside step function. In a lattice model with coordination number z the 
dwell time at a site with m bridges is calculated by averaging over the hopping 
rates to the nearest neighbors (see Fig. 3.4). 



oo oo ^ 

Pn.,---Pm.^^, ; (3.10) 



Fortunately, this ensemble averaging procedure can be accurately approximated by 
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an effective Arrfienius relation: 

= Toexp[A(m - m)] (3.11) 

Tfie validity of the approximation is most important for sites with m>m bridges, 
since the diffusive exploration allows the particle to quickly cascade into these deep 
energy wells. 

The effective Arrhenius relation greatly simphfies the calculation, since so long 
as Am is sufficiently large, the probability of the particle still being attached to the 
surface after an n step random walk is (1 — X^r^)"'"^ — \^~ exp(— Am)]"~^. Here 
Km — ^ exp(— Am) is the departure rate from a site with m bridges. Interestingly, 
in this approximation scheme the attachment probability is independent of the 
particular bridge numbers {mi,...,m„} realized during the walk. Thus, the 
probability of departure /„ after exactly n steps is: 

/„ = [exp(7)-l]exp(-7n) (3.12) 
7 = - log[l - exp(-Am)] (3.13) 

The average number of steps for the random walk is J^^i^/n = exp(Am). To 
calculate the departure time distribution $(i) we use /„ to average over the departure 
time distribution for walks with a given n, 0„(t). 

oo 

m = Y.f-^n{t) (3.14) 

n=l 

«o^n(|/.f-.(^))^(-|:'*) (-) 

Here Sjn{t) is the survival probability at time t for a site with m bridges, used to 
determine the probability of departure between t and t + dt. If there was only one 
hopping pathway with rate k, we would have — K,exp{—Kt). The generalization 
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accounts for the fact that the particle can hop to any of its z neighbours, and the 
probabihty of departure is not simply exponential. 



'S'm(^) =wliPa exp {-tK^-^c\ 



(3.16) 



va=l 



It is convenient to Fourier transform so that one can sum the resulting 

geometric series for $(a;). 



/oo 
0„(t)expHu;i]di = X(a;)' 
-oo 



X{u) = ^ ^ • • • ^ • • -P, 



m=l mi=l mj=l 



i=l 



i=l 



^u) = [exp(7) - 1] 5^ [exp(-7)X(u;)]" = [exp(7) - 1] 



X(uj) 



n=l 



exp (7) — X{uj) 



(3.17) 
(3.18) 
(3.19) 



To facilitate a simpler calculation, we employ a coarse-graining procedure to 
dispense with the tensor indices {m, mi, ...,171^} in the definition of X{u!). In the 

z 

summation there are many terms for which the value of ^ K^-^mi a-re equal, but with 

_ _ _ i=i ' 

different weight factors PmPmi ■ ■ ■ Pm^- To eliminate this degeneracy we introduce a 
smooth function /(«;) normalized according to J f{K)dK = 1. 



X{u;) / f{K 



K 

I , 

K-\-iuj 



-dK 



(3.20) 



The inverse Fourier transform is performed using the residue theorem to obtain 
the final result. The contour integral is closed in the upper half plane, with all the 
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poles on the imaginary axis ai u — iz. 



27r 



27r 



reSu 



r=\ 



exp[ici;t]X(u;) 



= [exp(7) - Iji^res^ 

r=l 

oo 

= [exp(7) - l]i 



r=l 



exp(7) — X{uj) 

(3.21) 

[lujrt] [XM + (u; - g;.) (f ),^,^ + • • • } 
exp(7) - {XK) + (a; - c,) (f + • • • } 

e'X.\)[iuJrt\X{(jJr 



exp 



= exp(7)[exp(7) - 1] ^ 



exp(— z^t) 



ZrY 



(3.22) 



Here 2;^ labels the roots of the equation 



exp(7) — X{iz) — 



(3.23) 



The benefit of the coarse-graining is now more transparent, as the residues are all 
labeled by a single index r as opposed to the tensor indices {m,mi, m^}. 

In Fig. 3.6 the departure time distribution is plotted in the diffusive regime. 
The optimal regime for fast particle departure is to have a large number (m ~ 10) of 
weakly bound key-lock bridges. In this scenario the departure time distribution is 
accurately approximated as a single exponential, $(t) = Kjnexp{—Krnt). 

We now discuss the behavior of the departure time distribution in several regimes 
of interest. At fixed m, for small A the behavior is non-universal. The departure 
time distribution exhibits multi-stage behavior, where the initial departure and 
long time behavior may both take the shape of a power law, albeit with different 
exponents (see A = 0.5, 1 curves in Fig. 3.5). 

As the strength of the key-lock binding increases (A > 1) there is a crossover 
from non-universal behavior to universal power law behavior for the first several 
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decades in time (see A = 2, 3 curves in Fig. 3.5). 

<l>(t) ~ t-'^ (3.24) 

For A > 3 we enter the regime of multiexponential beating. The initial 
departure behavior is well described as an exponential with initial departure rate 
= exp(— Am)/ro and characteristic timescale ~ 15ro. 

$(t) ~ Krnexp{-K*t) (3.25) 

We attribute n* to the diffusive cascade of particles from states with m bridges 
into more highly connected states. Since this process involves particles finding a 
lower energy state, k* does not depend on A. As indicated by the small departure 
probability, the binding is nearly irreversible in this regime (see A = 5, 7 curves in 
Fig. 3.5). 




Figure 3.5: Departure time distribution function versus time with m = 3. The dotted hues for 
A = 5 and 7 show the exponential approximation = ii'm exp(— 

We also plot the departure time distribution relevant to the experimental 
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Figure 3.6: Departure time distribution function versus time as determined by Eq. 3.21 in the 
diffusive regime. In the plot the average binding energy is held constant at 4T. The theoretically 
determined departure time distribution can be compared to an experiment with DNA-grafted 
colloids which observed power law behavior with exponent —1.5. 

situation where the average binding energy is held constant [30]. 



The optimal regime for fast departure is to have a large number (m ~ 10) of weakly 
bound bridges (see Fig. 3.6). In this fast departure regime the departure time 
distribution is well approximated as a single exponential, $(t) = Kj^ex^i^—Kmi) . 



We now turn to discuss the statistics for the in-plane diffusion of the particle. 
We first note that the in-plane trajectory of the particle subjected to a delta- 
correlated random potential remains statistically equivalent to an unbiased random 
walk. As a result, the mean squared displacement for an n step random walk 
remains (r^) = na^. As the particle explores the landscape it cascades into deeper 




1 — exp(— m) 



(3.26) 



3.6 Diffusion 
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Figure 3.7: Root mean squared displacement vs. time with Am = 4. The curves are calculated 
from the parametric equations 3.29, 3.30. 

energy wells, the hopping time increases, and the diffusion gets slower. In the limit 
n — >■ cx) the average hopping time can be determined from the equilibrium canonical 
distribution. For Poisson distributed bridge numbers m, this corresponds to a finite 
renormalization of the diffusion coefficient D* with Dq = a? /Atq. 



~ exp (me^) — 1 

t) = n {Trn) = UTo CXp (-Am) } P„ exp (Am) = UTq , . j—- -r 

exp (Am) [exp(m) — IJ 

(3.27) 

_ 1 g (r^) _ ^ exp (Am) [exp(m) - 1] 
^ =4^-^° exp (meA) - 1 ^^"^^^ 

This "ergodic" behavior is only achieved after a very long time. Generally, an 

n step random walk cannot visit sites with arbitrarily large m. In this transient 

regime one should only average over sites with m < m*. In the language of the 

statistics of extreme events, m* — 1 is the maximum "expected" value of m in a 

sample of n independent realizations [79]. Even with this complication, the average 

diffusion time (t) and the mean squared displacement (r^) can both be expressed in 
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Figure 3.8: The dimensionless diffusion coefficient D = jjy^^jj^ plotted against time. 



terms of m*, which defines their relationship in parametric form. 



1 - 



P(m, m*) 



(3.29) 



(3.30) 



D* \ 1 — exp[—me^ 
Here P{x, m*) = 7(0:, m*)/T{m*) = exp(— x) x'^/kl is the regularized lower 

incomplete F function. In the limit m* —>■ 00 we recover the renormalized diffusion 
relation (t) = (r^) /D*, although this occurs at very long, often unrealistic times. In 
the transient regime we expect anomalous, subdiffusive behavior. As indicated in 
Fig. 3.7, this subdiffusive behavior is typical for strong enough key- lock interactions. 
Figure 3.8 is a plot of the dimensionless diffusion coefficient versus time. 

By approximating the incomplete gamma functions, the transient behavior may 
be well described by a power law with a single free parameter /3 ^ 0.15 (see Fig. 
3.9). 

m" (3.31) 

^0/ 



2xh 
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(3.32) 



A 

Figure 3.9: Comparison of the power law exponent determined numerically to Eq. 3.32 in the 
transient regime. 

There is an analogy between our results and dispersive transport in amorphous 
materials. In these systems a length dependence of the effective mobility [82] can 
be interpreted within the context of the statistics of extreme events [79]. The time 
required for charge carriers to travel through the material depends on the dwell times 
spent at all of the trapping centers. Since this transit time will be dominated by the 
dwell time of the deepest trapping center, one would like to know how the thickness 
of the material effects the distribution for the largest trapping depth. The analogy 
to our result is made by replacing the material thickness with the number of steps 
in the random walk n, and replacing the distribution for the largest trap depth with 
the distribution for m*, since the bridge number is related to energy by Eq. 3.1. 
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Figure 3.10: Graphical depiction of key-lock binding between nanoparticles functionalized with 
complementary ssDNA. The resulting structures can be disordered, fractal-like aggregates, or 
crystalline. 

3.7 Connection to Experiments 

We now would like to make a connection between our results and recent 
experiments with DNA-grafted colloids. The departure time distribution can 
be compared to an experiment which determined the time-varying separation of 
two DNA-grafted colloids in an optical trap [26]. In the experimental setup, two 
particles are bound by DNA bridges, and after breaking all connections diffuse to 
the width of the optical trap. Because the length of the DNA chains grafted on 
the particle is much shorter than the particle radius, surface curvature effects can 
be neglected. The interaction resembles that of a particle interacting with a patch 
on a 2D substrate. Experimentally the tail of the departure time distribution was 
observed to be a power law $(t) ~ t^^'^. Qualitatively similar behavior is predicted 
by the theory with m ~ 1 and average binding free energy of several T (see m = 1 
curve in Fig. 3.6). 
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Figure 3.11: Plot of the characteristic times Tdep and Tdif versus m at constant binding energy 
{const = 4 in Eq. 3.34). 

In addition, our work provides insight into the slow crystalhzation dynamics of 
key-lock binding particles (see Fig. 3.10). In [1], l/im diameter particles grafted 
with ssDNA formed reversible, disordered aggregates. The average number of 
key-lock bridges between particles was m ~ 2. The authors of [26] observed random 
hexagonal close packed crystals by further reducing the surface density of DNA 
strands on the particles. The crystallization process requires that particles rapidly 
detach and reattach at the desired lattice location. In the localized regime particle 
desorption is the relevant process. 

In the diffusive regime surface diffusion also plays a role in the rearrangement of 
particles into the desired crystalline structure. To determine which process is more 
important for particle rearrangement, we can compare the departure time with the 
time required for a particle to diffusively explore the surface of a particle to which it 
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is bound. The time T^ep required for 90% of the particles to depart is: 




(3.33) 



To estimate the time required for diffusive rearrangement T^if we use the 
parametric equations 3.29, 3.30. In [1] particles of radius R = were grafted 

with DNA chains of length I ~ 20nm. Assuming the correlation length a ~ Z we 
have ~ (^)^ — 10^- Figure 3.11 shows a comparison of Tdif and Tdep at 
constant binding free energy. 

A777 

7 — + log(l — exp(— mj) = const (3.34) 

1 — exp(— m) 

This expression for the binding energy takes into account the entropy reduction 
associated with the non-ergodic degrees of freedom. For a detailed discussion of this 
topic see reference [30]. Since T^if > T^ep, colloidal desorption and reattachment is 
the dominant mechanism by which particles rearrange. 

As the figure indicates, the optimal regime of fast departure is to have a large 
number (m > 10) of weakly bound key-lock bridges. We predict a localized regime 
below the crossover where particle departure is relatively fast. Just beyond the 
crossover there is a relative maximum in T^ep before it decreases at large m. The 
increase in departure time at the onset of diffusive behavior is indicative of a regime 
where the system ages. In this regime the interplay of diffusion an desorption leads 
to longer bound state lifetimes, and an increase in the departure time. 

We now turn to the question of designing a future experiment which will facilitate 
fast particle departure and hence colloidal crystallization. The essential goal is to 
increase the average number of key-lock bridges between particle pairs. Increasing 
the surface density of DNA strands alone results in the formation of a brush, which 
decreases the effective cross section for the interaction between complementary 
DNA. Instead we propose the introduction of long, flexible DNA linkers between 
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Figure 3.12: Graphical depiction comparing recent experiments with DNA-grafted coUoids to a 
future implementation with long flexible linkers. Increasing the number of key-lock bridges 
between particle pairs potentially decreases the time required for crystallization. 

particles with a high coverage of short ssDNA (see Fig. 3.12). This system has 
the potential to realize more key-lock bridges between particle pairs as compared 
to previous experiments, and therefore substantially reduce the time required for 
crystallization. 

3.8 Summary 

In this chapter we studied the dynamics of particles which interact through 
the reversible formation of multiple key-lock bridges. Well before the percolation 
threshold is reached there is a crossover from a localized regime to a diffusive regime. 
In the localized regime the particles remain close to their original attachment site 
until departing. In this regime particles are attached to finite clusters, and the 
system exhibits an exponential distribution of departure times. Once the radius of 
gyration of the cluster exceeds the characteristic radius for the particles' random 
walk, the finite clusters behave effectively as infinite clusters. Diffusion allows the 



63 



particles to cascade into deeper energy wells, which leads to a decrease in hopping 
rate. The diffusion slows and the bound state lifetime increases, a phenomenon 
qualitatively similar to aging in glassy systems. In the diffusive regime we discussed 
the statistics for the particles' in-plane diffusion. Weak key-lock interactions 
give rise to a finite renormalization of the diffusion coefficient. However, as the 
strength of the interaction increases (larger A), the system exhibits anomalous, 
subdiffusive behavior. This situation is analogous to dispersive transport in 
disordered semiconductors. We then made the connection between our calculation 
of the departure time distribution and recent experiments with DNA-coated colloids. 
The findings indicate that the optimal regime for colloidal crystaUization is to have a 
large number of weakly bound key-lock bridges. A modified experimental setup was 
proposed which has the potential to realize this regime of fast particle departure. 
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CHAPTER IV 

SELF-ASSEMBLY OF DNA-CODED CLUSTERS 
4.1 Motivation and Problem Description 

Over the past decade, a number of proposals have identified potential 
applications of DNA for self-assembly of micro- and nanostructures [18], [17], 
[46], [47], [48]. Among these proposals, one common theme is finding a way to 
utilize the high degree of selectivity present in DNA-mediated interactions. An 
exciting and potentially promising application of these ideas is to use DNA-mediated 
interactions to programmable self-assemble nanoparticle structures [25], [15], [45], 
[44]. Generically, these schemes utilize colloidal particles functionalized with 
specially designed ssDNA (markers), whose sequence defines the particle type. 
Selective, type-dependent interactions can then be introduced either by making the 
markers complementary to each other, or by using linker-DNA chains whose ends 
are complementary to particular maker sequences. Independent of these studies, 
there are numerous proposals to make sophisticated nano-blocks which can be used 
for hierarchical self-assembly. One recent advance in the self-assembly of anisotropic 
clusters is the work of Manoharan et. al [3]. They devised a scheme to produce 
stable clusters of n polystyrene microspheres. The clusters were assembled in a 
colloidal system consisting of evaporating oil droplets suspended in water, with the 
microspheres attached to the droplet interface. The resulting clusters, unique for 
each n, are optimal in the sense that they minimize the second moment of the mass 
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distribution M2 = X]r=i('^* " "^cmf- 

In this chapter [61], we present a theoretical discussion of a method which 
essentially merges the two approaches. We propose to utilize DNA to self-assemble 
colloidal clusters, somewhat similar to those in Ref . [3] . An important new aspect of 
the scheme is that the clusters are "decorated": each particle in the resulting cluster 
is distinguished by a unique DNA marker sequence. As a result, the clusters have 
additional degrees of freedom associated with particle permutation, and potentially 
may have more selective and sophisticated inter-cluster interactions essential for 
hierarchic self-assembly. In addition, the formation of such clusters would be an 
important step towards programmable self-assembly of micro- and nanostructures of 
an arbitrary shape, as suggested in Ref. [65]. 




n=7 n=9 

Pentagonal dipyramid Tnaugmented tnangular prism 




Figure 4.1: The minimal second moment clusters for rt = 5, 6, 7, and 9. Pictures of all the clusters 
from ?T. = 4 to 15 are available in [3]. 

We begin with octopus-like particles functionalized with dsDNA, with each 
strand terminated by a short ssDNA marker sequence. We assume that each particle 
i has a unique code, i.e. the maker sequence Sj of ssDNA attached to it. We then 
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Table 4.1: Minimal Second Moment Clusters 


n 


Polyhedra Name 


Schonflies Point Group 


4 


Tetrahedron 


a 


5 


Triangular dipyramid 




6 


Octahedron 


Oh 


7 


Pentagonal dipyramid 




8 


Snub disphenoid 




9 


Triaugmented triangular prism 




10 


Gyroelongated square dipyramid 


Dm 


11 


(non-convex) 


Cs 



introduce anchor DNA to the system, ssDNA with sequence SASB---Sn, with Sj the 
sequence complementary to the marker sequence Si. The anchor is designed to 
hybridize with one particle of each type. Consider a cluster of n particles attached 
to a single anchor. 




Figure 4.2: Schematic representation of the method for constructing decorated colloidal clusters 
using ssDNA " anchors" . 

If we treat the DNA which link the particles to the anchor as Gaussian chains, 
there is an entropic contribution to the cluster free energy which can be expressed in 
terms of the particle configuration {ri, ...,r„} as follows. Here Rg is the radius of 
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gyration of the octopus-like DNA arms. 



3T 
2^ 



n 



F 



EC 




(4.1) 



'9 i=l 



This approximation of the DNA arms as Gaussian chains is acceptable provided 
their length L exceeds the persistence length 1^ ~ 50 nm and the probability of 
self-crossing is small [55]. The physical mechanism which determines the final 
particle configuration in our system is quite different from the capillary forces of 
Manoharan et al. However, because the functional form of the free energy is 
equivalent to the second moment of the mass distribution, the ground state of the 
cluster should correspond to the same optimal configuration. 



Consider a system with n particle species and an anchor of type 'saSb-'-'Su- 
The clusters we would like to build contain n distinct particles (each particle in the 
cluster carries a different DNA marker sequence) attached to a single anchor. Let 
Cn denote the molar concentration of the desired one anchor cluster. Because there 
are many DNA attached to each particle, multiple anchor structures can also form. 
The question is whether the experiment can be performed in a regime where the 
desired one anchor structure dominates, avoiding gelation. 

We consider the stability of type C„ with respect to alternative two anchor 
structures. To do so we determine the concentration Cn+i of 1 particle structures 
which are maximally connected, but do not have a 1 : 1 : • • • : 1 composition. In 
particular, these structures contain more than one particle of each type, which could 
cause problems in our self-assembly scheme [65] . There are also n particle structures 
Cn with the correct composition, but which contain two anchors. We would like to 
avoid the formation of these structures as well, as their presence decreases the overall 
yield of type C„. Figure 4.3 enumerates the various structures for an n = 3 species 



4.2 Equilibrium Treatment 
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system. If the experiment can be performed as hoped, we will find a regime where 



the ratios 



c, 



n + l 



Cn 



and are small. To this end, the equilibrium concentrations C„, 



Cn+i, and Cn are determined by equilibrating the chemical potential of the clusters 
with their constituents. 




particle of type k 



anchor 



I DNA linker 



Figure 4.3: The topologically distinct one and two anchor structures for an anchor ssDNA with 
sequence 'sa'sbsc- Different structure varieites may be obtained by relabehng the particle indices 
subject to the constraint that no more than one particle of each type is attached to a given anchor. 



First we determine the molar concentration C„ of the desired one anchor 
structure, composed of one anchor and n particles. Let q denote the molar 
concentration of species and Ca the molar anchor concentration. Here Co = IM is 
a reference concentration. 



Tlog 
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i=l 
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i=l 
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(4.2) 

(4.3) 
(4.4) 



Here e„ is the binding free energy of the cluster, which has a contribution from 
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the hybridization free energy AGj associated with attaching particle i to the anchor, 
and an entropic contribution from the number of ways to construct the cluster (since 
each particle has N hybridizable DNA arms). In addition we must take into account 
the entropy for the internal degrees of freedom in the structure stemming from the 
flexibility of the DNA attachments to the anchor. In the Gaussian approximation, 
neglecting the excluded volume between particles this localization volume v '-^ Rg 
can be calculated exactly. 

3T 



F(ri,...,r„) = 7^5]^' 
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d?ri...d^rn exp 



T 



9 \ 3 



(4.5) 
(4.6) 
(4.7) 



We now consider the competing two anchor structure Cn+i- To determine the 
equilibrium concentration it is instructive to consider the reaction in which two 
clusters of type C„ combine to form a single cluster of type C„+i and release n — 1 
particles into solution. Since there are many DNA attached to each particle, in 



what follows we omit factors of "^^y^- 
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(4.8) 
(4.9) 
(4.10) 



i=l 



The localization volume V2 can be calculated in a similar fashion, fixing one particle 
at the origin and integrating over the n remaining particle positions {ri, ...,r„} and 
the position of the two anchors and ra2- We make use of the following formula 
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for multivariate Gaussian integrals. 
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F(ri, r„i, r„2) = ^[r^, + (r„ - r^a)' + E {(r. - r^i)^ + (r„2 - r,)^}] (4.12) 
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Similarly one can obtain the cluster concentration for the two anchor structure C„ 
The localization volume ^3 can be calculated in a similar fashion to v^- 



V3 = Rt 
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(4.15) 



We consider the symmetrical case AGj = AG and equal initial particle 
concentrations cf^ — c^°^ for all species i. In this case we have CA — CB — --- — Cn = c. 
We can express the concentration of the competing two anchor structures Cn+i and 
Cn in terms of the concentration of the desired one anchor cluster C„. Since there 
are many DNA attached to each particle, we omit factors of 

AG' 



Cn = Ca yNvc exp 

Cn+1 ^ 



2 



V^" C" 



(4.16) 

(4.17) 

(4.18) 



v^" c" 

The concentration of free anchors can be determined from the equation for anchor 
conservation. 

= Ca + Cn + 2Cn + 2nC„+i (4.19) 
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Figure 4.4: The molar concentrations Cq, C„, nC„+i, and C„ in tire symmetrical case for a system 
with n = 5 particle species. The total particle volume fraction ncj) w .25 and -r^ — 10 . 



We are interested in the low temperature regime where there are no free anchors 
in solution. We determine the saturation values for the ratios of interest by noting 
that the Boltzmann factor 5 = exp [— ^] ^ 1 in this regime. 



(4.20) 



(4.21) 



Since 



Cn-t 



^ 1, eq. 4.21 provides the experimental constraint for suppressing 



the two anchor structures. Taking the radius of the hard spheres R ^ Rg, it can 
be interpreted as a criterion for choosing the initial anchor concentration ci"'' for an 



n 



species system with = ^R'^c^°^ the particle volume fraction for an individual 



species. 
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(4.22) 



72 



The condition gives the maximum anchor concentration for the two anchor structures 

(o) 

to be suppressed. Since < ^ the theoretical hmits are ^ < 1, .29, .06, and 
.01 for n — 4, 5, 6, and 7 respectively. In figure 4.4 we plot the solution for the 
concentrations. There is a large temperature regime < 2) where the two anchor 
structures are suppressed in favor of the desired one anchor structures. 

4.3 Irreversible Binding 

In the previous section we performed an equilibrium calculation to determine 
the yield of the cluster C„ [61]. The results of that study indicated that the 
concentration of anchors must be kept very small to prevent the aggregation of larger 
clusters (see Eq. 4.22). From an experimental perspective this result is somewhat 
disappointing, since the overall yield of the cluster C„ is proportional to the anchor 
concentration. The situation is considerably improved in the regime of irreversible 
binding of particles to anchors. In what follows we present a calculation for the yield 
of the desired one anchor cluster far from equilibrium. To distinguish between the 
results of the previous section and the regime of irreversible binding we will change 
the notation slightly. This change reflects the fact that the role of the DNA anchor 
could also be played by a patchy colloidal particle as discussed below. Henceforth 
we will refer to the DNA anchors as DNA scaffolds. In the new notation the desired 
one anchor cluster C„ is called the star cluster. 

The plan for this section is the following. The goal is to maximize the yield 
for the star cluster. We analytically calculate the yield of the star cluster in the 
regime of irreversible binding. The analytical results are compared to the numerical 
results for the full aggregation equations. Prom an experimental perspective, the 
most important result is the determination of an optimal concentration ratio for 
experiments (see Eq. 4.30). To conclude we discuss the experimental feasibility of 
the self-assembly proposal. 
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The basic idea behind the procedure is as follows (see Fig. 4.5). Particles are 
functionalized with single-stranded DNA (ssDNA) markers which determine the 
particle color. There may be many DNA attached to each particle, but on any 
given particle the marker sequence is identical. One then introduces DNA scaffolds 
to the system. The scaffold is a structure with n ssDNA markers, each marker 
complementary to one of the particle colors. Hybridization of the ssDNA markers 
on the particles to those on the scaffold results in the formation of colored particle 
clusters. Because there are many DNA attached to each particle, clusters can form 
which contain more than one scaffold. The essential goal of the procedure is to 
maximize the concentration of a particular type of cluster which we denote the star 
cluster. The star cluster contains one and only one scaffold to which n particles are 
attached, each particle having a distinct color. 

We should note that the role of the scaffold could also be played by a patchy 
particle ([83], [84], [85]). For example, these patches are regions on the particle 
surface where one can graft ssDNA markers. In this case there may be several DNA 
connections between a patch and colored particle. Our conclusions will still be 
valid, provided the patch size is chosen so that a patch interacts with at most one 
particle. 

To understand the basic physics behind the aggregation process we consider the 
mobility mismatch between the particles and the scaffolds. In solution, a particle 
with radius R ~ l/im has a diffusion coefficient given by the Stokes- Einstein relation 
D — kBT/dnrjR. On the other hand, the size of the scaffold a ~ lOnm. As a 
result the scaffolds diffuse R/a ~ 100 times faster than the particles. To first 
approximation the resulting aggregation is a two stage process. In the first stage 
the particles recruit different numbers of scaffolds via the fast scaffold diffusion and 
subsequent DNA hybridization. Since we consider the regime of strong binding 
where these bonds are irreversible, the result is a Poisson distribution over the 



74 







DNA scaffold 



Colored particle 



Figure 4.5: A graphical depiction of the scheme for self-assembhng star clusters using DNA 
scaffolds. In the diagram (not drawn to scale) the scaffold funcionality 71 = 4. 

concentration of particles with m scaffolds attached. Let Ci denote the concentration 
of particles with color i, and c denote the total concentration of scaffolds. The total 
particle concentration Ctot = Y^^=i ^i- The concentration C^"^^ of particles of color i 
with m scaffolds attached is 



P 



a 



p"^ exp(— p) 



ml 



(4.23) 
(4.24) 



tot 



In the second stage there are no free scaffolds left in solution, and these particles 
decorated with scaffolds aggregate to form the final clusters. The seed to build a 
star cluster is a particle of any color with exactly one scaffold attached. This seed 
must aggregate with n — 1 particles of different colors, each of which has no scaffolds. 
We now calculate the concentration of the star cluster C*. The yield of the desired 
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star cluster is quantified in terms of the star mass fraction = {nC^) /Ctot- 



Here x — np is the scaffold functionality n multiplied by the concentration ratio 
p. By choosing p — 1/n the mass fraction attains a maximum of exp(— 1) ~ 0.37! 
This result indicates that by selecting the appropriate scaffold concentration, in the 
nonequilibrium regime up to 37% of the particles will aggregate to form star clusters. 
This is a significant improvement over the situation in the equilibrium regime. 

This treatment of the problem captures the physics of star cluster formation, but 
it does not account for the loss of star clusters due to aggregation. In particular, 
as long as there are scaffolds with markers available for hybridization, when these 
scaffolds encounter a star cluster they can aggregate to form a larger cluster. We 
now estimate how this aggregation effects the final concentration of star clusters. 

Consider the beginning of the second stage in our aggregation process. There 
are no longer any free scaffolds in solution, but a scaffold can have up to n — 1 
DNA markers still available for hybridization. We would like to determine how the 
star cluster mass fraction M^{y) changes as a function of the fraction of saturated 
scaffolds y. Here a saturated scaffold has particles hybridized to all n of its DNA 
markers, and is therefore unreactive. If s is the expectation that a slot on the 
scaffold is filled, then the fraction of saturated scaffolds \s y — The average 

number of open slots on a scaffold is (n — 1)(1 — s). Consider filling an open slot 
on the scaffold. The probability that the particle which filled the slot was part of 
a star cluster is M*(y). The average rate r{y) at which star clusters are lost to 
aggregation is then 



Here the exponent a = {n — 2)/{n — 1). We can then construct a differential 




(4.25) 
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(4.26) 
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equation for taking into account this loss due to aggregation. 



xr(y) (4.27) 



dy dy 

In the absence of this loss term the result of the calculation should recover our previous 
result Eq. 4.25. This zeroth order approximation is just Mi°\y) = xy exp(—xy) 
which gives the correct star cluster concentration once all of the scaffolds are 
saturated {y — 1). To simplify the analysis a bit we take a — 1 which is an 
excellent approximation in the limit of large scaffold functionality n. This is an 
inhomogeneous first order differential equation which can be solved by introducing 
an integrating factor u{y) — y^. The initial condition which must be satisfied is 
M*(0) = 0. We are interested in the final star mass fraction M*, which is M^{y = 1). 
The result is 
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(4.28) 



Here r(x) is the gamma function and E^{x) = ft exp{—xt)dt is the exponential 

1 

integral of order u. 

We can perform a similar type of analysis in the case when there is only one 
particle color. In this case the n ssDNA markers on the scaffold all have identical 
sequences complementary to this color. It turns out that the result for the mass 
fraction is the same. Because the mass fraction is the same in both cases, we can 
gain insight into the behavior of the system with many colors by analyzing the much 
simpler one color system. To test our predictions, we numerically solved a system 
of differential equations which models the irreversible aggregation between particles 
(one color) and scaffolds. 

= - ^ KjjjiyCijCj'ji - Cij KjjijCjj (4.29) 

i+i'=I i,j 
j+j'=J 
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This equation is the Smoluchowski coagulation equation [86] adapted to our 
system. Cij is the concentration of the cluster with i scaffolds and j particles. 
Kijiiji is the rate constant for the irreversible reaction Cij + Cj/j/ — > Cj+j/j+j'. We 
assume that the rates are diffusion limited in which case we can estimate the rate 
for any pair of clusters by Kijiij/ = iirDgRi. The larger cluster with hydrodynamic 
radius Ri ~ n^^ plays the role of a sink. Here rii is the number of particles in the 
larger cluster and Dg — 1 6nr]Rs is the diffusion constant for the smaller cluster. 
To simply matters we only consider tree like structures, i.e. we do not consider the 
formation of clusters with internal loops. We have truncated the set of equations by 
considering clusters with a maximum of 10 scaffolds. 

By solving these equations we can determine the concentration of stars = Ci„ 
in this notation and test the validity of our two stage ansatz. As indicated in Fig. 
4.6, the result of our analytical calculation matches the results of the full numerical 
calculation up to an overall normalization factor of order unity. Several points are 
in order. 

The optimal concentration ratio p for experiments is easily determined from 
= 0. The result is Xmux — 0.47. For scaffolds of functionality n the 
concentration ratio should be chosen as: 

p=^. (4.30) 

n 

Note that the maximum attainable star cluster yield M*(a;niax) 1/4 does not 
decrease with increasing n. In fact, the numerical results predict a slight increase 
in star cluster yield for larger n. Solving the aggregation equations becomes 
computationally expensive, but it can still be done by reducing the maximum 
number of scaffolds in a cluster. For example, considering clusters with up to 
5 scaffolds for n = 10 gives M*(a;inax) — 0.3. These results are important from 
the perspective of experimental feasibility for the self-assembly method. This is 
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Figure 4.6: The mass fraction as a function of p for scaffolds with functionahty n = 3 (red), 4 
(blue), 5 (green), and 6 (black). The results determined numerically from the full solution of the 
Smoluchowski coagulation equation (markers) can be compared to results of the anlaytical 
calculation (lines) Eq. 4.28. The resulting agreement is good up to an overall normalization factor 
7„ in the range 1.2 to 1.5 which normalizes the analytical curves. 

to be contrasted with the earher equihbrium treatment. There the condition to 
suppress the aggregation of larger clusters imposed a fairly strict constraint [61] on 
the concentration ratio p < n^^^ («)" ^- From the perspective of self-assembling 
stars with large n this renders the regime of irreversible binding far more appealing 
than the equilibrium regime. 

If an experiment is performed with the optimal concentration ratio, the clusters 
which self-assemble are easily separated by density gradient centrifugation[87]. In 
this regime most of the particles are monomers, in star clusters, or in saturated two 
scaffold clusters. These clusters contain, 1, n, and 2n — 1 particles respectively. 
The disparity in hydrodynamic radius and sedimentation velocity of these clusters 
makes the separation procedure experimentally feasible. 

In this section we considered a DNA scaffold method for self-assembling star 
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clusters of n colored particles. By taking advantage of the mobility mismatch 
between particles and scaffolds, we were able to formulate a nonequilibrium 
calculation of the star mass fraction. The results of the calculation were compared 
to the numerical results of the full Smoluchowski coagulation equation for the system. 
Good agreement is established between the analytical calculation and the numerics. 
In the regime of irreversible binding the yield of the desired star cluster is drastically 
improved in comparison to earlier equilibrium estimates. In nonequilibrium we 
find an experimentally feasible regime for the self-assembly of star clusters with a 
maximum mass fraction ~ 1/4. We determined the optimal concentration ratio 
for an experimental implementation of our proposal. The additional color degrees 
of freedom associated with particle permutation in these clusters makes them ideal 
candidates as building blocks in a future hierarchical self-assembly scheme. In 
addition, these clusters can serve as the starting point to self-assemble structures of 
arbitrary geometry[65] . The experimental realization of self- assembling star clusters 
using DNA scaffolds would constitute an important step towards realizing the full 
potential of DNA mediated interactions in nanoscience. 

4.4 Cluster Degeneracy 

Building these decorated colloidal clusters is the first major experimental step in 
a new self-assembly proposal which will be discussed in the next chapter. However, 
in order to utilize the resulting clusters as building blocks, an additional ordering 
is necessary. The problem is that the decoration introduces degeneracy in the 
ground state configuration. This degeneracy was not present in [3] since all the 
polystyrene spheres were identical. Namely, in the colloidal clusters self-assembled 
by our method, permuting the particle labels in a cluster does not change the second 
moment of the mass distribution (see Figure 4.7). We need a method to select a 
single "isomer" out of the many present after self-assembly. In the DNA-coUoidal 
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system considered here, this isomer selection can be facihtated by "hnker" ssDNA. 
These are short ssDNA with sequence saSb to connect particles A and B. We first 
construct a list of nearest neighbors for the chosen isomer, and introduce linker DNA 
for each nearest neighbor pair. The octopus-like DNA arms of the given particles 
will hybridize to the linkers, resulting in a sping-like attraction between the selected 
particle pairs. Note that the length L of the DNA arms must be on the order of 
the linear dimension of the original cluster. Otherwise the interparticle links cannot 
form upon introduction of linker DNA to the system. It should be noted that 
although this method breaks the permutation degeneracy of a cluster, the right-left 
degeneracy will still be present. 





Figure 4.7: An illustration of degeneracy in DNA-coded nanoclusters. Two different n = 6 isomers 
are pictured, both with the same minimal second moment configuration, the octahedron. 
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CHAPTER V 



PROGRAMMABLE SELF-ASSEMBLY 
5.1 Introduction 

Over recent years, significant attention has been attracted to the possibihty 
of nanotechnological apphcations of DNA [18], [17], [46], [47], [48]. Among the 
various proposals, one of the most interesting directions is the use of DNA-mcdiatcd 
interactions for programmable self-assembly of nanoparticle structures [25], [15], [45], 
[44] . Several schemes of such self-assembly have been studied both experimentally 
and theoretically. Their common theme is the use of colloidal particles functionalized 
with specially designed ssDNA (markers), whose sequence defines the particle type. 
In such systems, selective type-dependent interactions can be introduced either by 
making the markers complementary to each other, or by using linker-DNA chains 
whose ends are complementary to particular marker sequences. 

Recent theoretical studies have addressed the expected phase behavior [56], 
melting properties [2], and morphological diversity [55] of DNA-coUoidal assemblies. 
In particular, there are indications that these techniques can be utilized for fabrication 
of photonic band gap materials [74], [75]. Despite significant experimental progress, 
the long-term potential of DNA-based self-assembly is far from being realized. For 
instance, most of the experimental studies of DNA-coUoidal systems report random 
aggregation of the particles [53] . Some degree of structural control in these systems 
has been achieved, mostly by varying the relative sizes of particles, rather than by 
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tuning the interactions [43], [16]. 

In the present chapter [65] , we take a broader view of programmable self-assembly. 
While this theoretical study is strongly motivated by the prospects of DNA-coUoidal 
systems, our main objective is to address a more general question: how well can a 
desired structure be encoded by tunable interactions between its constituents? The 
particular model system on which we focus consists of distinguishable particles with 
individually controlled interactions between any pair of them. In the first stage of 
our study, we analyze a simplified yet generic version of such a system. All the 
particles have the same repulsive potential, while the attraction is introduced only 
between selected pairs of particles in the form of a spring-like quadratic potential. 
This potential mimics the effect of a stretchable DNA molecule whose ends can 
selectively adsorb to the particular pair of particles. We then introduce a number of 
additional features which make the model a more realistic description of an actual 
DNA-coUoidal system. 

Our major result is that a combination of stretchable interparticle linkers 
(e.g. sufficiently long DNA), and a soft repulsive potential greatly reduces (or 
totally eliminates) the probability of self-assembling an undesired structure. The 
experimental prototype is a system of particles in a mixture of two types of DNA 
molecules which can selectively adsorb to the particle surface. The first type are 
DNA molecules with two "sticky" ends, i.e. both end sequences of the DNA are 
complementary to the particle marker sequence. With one end adsorbed to the 
particle surface, the remaining sticky end makes it possible to introduce an attractive 
interparticle potential between selected particle pairs. The second type arc DNA 
molecules with one sticky end which adsorbs to the particle surface. These DNA 
strands give rise to a soft repulsive potential of entropic origin between all particle 
pairs. 

There is a natural analogy between our problem and the folding of proteins 
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where interactions between amino acids encode the overall structure. However, 
it should be emphasized that the task of programmable self-assembly in colloidal 
systems is even more demanding than protein folding: any intermediate metastable 
configuration has a much longer lifetime and therefore means a misfolding event. 
Because of this, we were looking for a self-assembly scenario which does not 
require thermally activated escape from a metastable configuration. This makes the 
problem additionally interesting from the theoretical perspective of "jamming", a 
phenomenon actively studied in the context of granular and colloidal systems. Our 
results can be interpreted in terms of a jamming-unjamming transition controlled by 
the interaction parameters. 

It should be emphasized that the goal of this work is primarily conceptual, as 
opposed to providing a manual for the immediate experimental realization of ordered 
colloidal structures. Nevertheless, future experimental schemes will be forced to 
overcome obstacles presented by colloidal jamming. With this in mind, one of the 
most salient features of our model is the ability to smooth the energy landscape by 
tuning the interactions between particles. 

The plan for the chapter is as follows. In section 2, we address the problem within 
a simplified generic model which mimics the nanoparticle system with stretchable 
DNA connections. An unexpected and very encouraging result of this study is that 
the misfolding (or jamming) in the model system can be completely avoided for a 
certain set of parameters. In section 3 the original model is adapted to a more 
realistic situation which incorporates the random character of the DNA-mediated 
interactions. In section 4 we discuss the prospects for the future experimental 
implementation of our scheme. In section 5 we summarize the major results. 
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5.2 Beads and Springs Model 



Consider an isolated group of repulsive particles linked via a polymer spring to 
their desired nearest neighbours. We assume that the DNA marker sequence is the 
same for any two markers attached to the same particle, but different particles have 
different marker sequences (i.e. each particle has a unique code). In this case, the 
attraction between any two particles can be effectively switched on by adding DNA 
"linkers" whose ends are complementary to the corresponding marker sequences of 
the particles. As a first approach to the problem, we introduce a generic "Beads 
& Springs" model which incorporates essential features of the DNA-nanoparticle 
system. The model system contains N particles with pairwise (type-independent) 
repulsive potential U (r) . In general, this repulsion may have a hard-core or soft-core 
behavior, or be a combination of the two. In order to model the DNA-induced 
type-dependent attraction, we introduce a harmonic potential (linear springs) which 
acts only between selected pairs of particles [22]. Thus, the model Hamiltonian has 
the following form: 



Here a,/3 are the particle indices, are their current positions, and k is the spring 
constant. The connectivity matrix element Ja/s may be either 1 or 0, depending on 
whether the two particles are connected by a spring, or not. Our goal is to program 
the desired spatial configuration by choosing an appropriate connectivity matrix 
Ja/3- A natural construction is to put Jap — 1 for any pair of particles which must he 
nearest neighbors in the desired cluster, and not to connect the particles otherwise 
(i.e. put Jap = 0). This construction assures that the target configuration is the 
ground state of the system. 

Note that our problem is somewhat similar to that of heteropolymer folding. In 
that case, the selective interactions between monomers (e.g. amino acids in protein) 




(5.1) 
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are responsible for the coding of the spacial structure of a globule. Our major 
concern is whether the kinetics of the system will allow it to reach the ground state 
within a reasonable time. Unfortunately, since the Brownian motion of a typical 
nanoparticle is relatively slow (compared to molecular time scales), it is unrealistic 
to expect that our system will be able to find the target configuration by "hopping" 
between various metastable states, as in the case of protein folding. However, our 
case is different because the attractive force grows with distance, as opposed to the 
short-range nature of heteropolymer self-interactions. As we shall see below, this 
difference is essential, making it possible for the system to reach the ground state 
without stopping at any metastable configuration. 

We have performed a molecular dynamics simulation of the above model by 
numerically integrating its Langevin equation: 

b-H^ = -VaH + ri, (5.2) 

Here b is the particle mobility. The thermal noise has been artificially suppressed in 
this study (i.e. ?7 = ). In other words, we have assumed the worst case scenario: 
once the system is trapped in a local energy minimum, it stays there indefinitely. 
The equations of motion were solved numerically by a first order Runge-Kutta 
method. First, we studied a system of A*" = 49 distinguishable particles in 2D, 
whose native configuration was a 7 x 7 square cluster (see Figure 5.1). Their initial 
positions were random, and the connectivity matrix was constructed according to 
the above nearest-neighbor rule. 

First, we studied the case when U{r) is a hard-core potential. More precisely, 
the repulsive force was given by a semi-linear form: 

he (r) = -^^^ = ^o{d - r)e{d - r) (5.3) 

Here © is the unit step function, and d is the diameter of the hard sphere. The 
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Figure 5.1: Programmable self-assembly in 2D, studied within the Beads & Springs model. 49 
particles are all distinct and connected with linear springs to encode the desired configuration (7x7 
square). The jamming behavior is observed for the case of hard spheres (a). However, the 
assembly of the target structure can be achieved if the repulsion is sufficiently "soft" (b). 

parameter determines the strength of the hard-core repulsion, and it does not 
affect the resuhs, as long as kq ^ k. In our simulations, we found that the hard 
sphere system eventually stops in a configuration definitely different from the desired 
one, a behavior which is well known in the context of granular and colloidal systems 
as "jamming" [88]. Remarkably, the jamming can be avoided when the hard-core 
repulsion is replaced by a soft-core potential: 

Usc{r) = Uoexp{-r/X). (5.4) 

Here the decay length A is of the order of the equilibrium interparticle distance tq. 
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This indicates that the energy landscape can be made smooth by a combination of 
long-range selective attraction and soft-core repulsion. 

The result is surprising and remarkably robust. In particular, in order to expand 
our finding to the 3D case, we studied the self-assembly of particles into tetrahedra 
of various sizes (A^ = 10, 20, 35). This time, the hard core interaction potential was 
superimposed with a soft shell repulsion, which makes the model more relevant for 
an actual DNA-coUoidal system: 



After the system has fully relaxed, a geometric measure of the folding success is 
determined by comparing particle separations of the desired final state to those 
generated from a set of random initial conditions. Figure 5.2 shows the "jamming 
phase diagram" for these systems. To assign a point on the diagram to the correct 
folding regime, we required 100 consecutive successful folds. While this criterion 
can only give an upper bound on the jamming probability (which is approximately 
1%), an additional analysis gives strong evidence that the correct folding region of 
the diagram corresponds to zero probability of jamming. 




(5.5) 
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Figure 5.2: "Jamming phase diagram" obtained for programmable self-assembly of tetrahedral 
clusters, within the Beads & Springs model. The control parameters depend on the equilibrium 
interparticle distance Tq, the diameter of the hard sphere d, and the range of the soft-shell 
repulsion A. 

5.3 Self-Assembly in DNA Colloidal Systems 



As we have seen, the introduction of a soft-core repulsive potential Use is crucial 
to a successful self-assembly proposal. In a real system, this repulsion can be 
generated e.g. by DNA or another water soluble polymer adsorbed to the particle 
surface. The mechanism is quite independent of the monomer chemistry, but for 
the sake of concreteness we will speak of the repulsion generated by DNA. Namely, 
we assume that a certain fraction of the DNA "arms" of the "octopus-like" particles 
are not terminated by a sticky end, and only play the role of a repulsive "buffer". 
When the polymer coverage is sufficiently low, the interparticle repulsion is primarily 
due to entropy loss of a chain squeezed between two particles. The characteristic 
length scale of this interaction is given by the radius of gyration of the "buffer" 
chain, Rg. The corresponding repulsive force fsc can be calculated exactly in the 
limit of relatively short buffer chains, Rg <^ d. The result of this calculation [55] can 
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be adequately expressed in the following compact form: 



fsc (r) 



ANRgT 
d{r-d) 



exp 



(r - d f 



(5.6) 



Here is the total number of buffer chains per particle. Note that even though 
this result is only valid for Rg <^ d, it correctly captures the Gaussian decay of the 
repulsive force expected for longer chains as well. Therefore, we expect the results 
to be at least qualitatively correct beyond the regime of short buffer chains. 

In addition to the modified soft potential, we have to take into account the 
random character of the realistic DNA-mediated attraction. It originates from the 
fact that (1) the number of DNA "arms" of the original octopus-hke particles will 
typically be determined by a random adsorption process, and (2) the fraction of 
the DNA chains recruited for linking a particular pair of particles is also random. 
In terms of our original model, this means that the "springs" will not have the 
same spring constant. If the individual linkers are modelled by Gaussian chains 
[20], the overall spring constant for a particular pair of connected particles is given 
by Kap — Tniap/^R'g , where Rg is the radius of gyration of a single linker, and 
map is the number of individual chains connecting the particles. We assume that 
this number obeys the generic Poisson distribution: P{m) — m"^e~"*/m!. As 
formulated, the model is cast as a system of coupled differential equations: 



Here Tap = \ra - rpl, h^p = (r« - r^) /r^/?. 

We have studied the behavior of the system as a function of two dimensionless 
parameters, one of which is the ratio of the buffer radius of gyration to the particle 
diameter, Rg/d . The other parameter characterizes the relative strength of the 




(5.7) 
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attractive and repulsive forces: 



m / R, 



a 



N \Rl 



(5i 



Including the above modifications of the model, the essential result of the study is 
that the jamming probability can be drastically suppressed, similarly to the previous 
case. However, the jamming in this system cannot be eliminated completely. 
Instead of an actual "phase boundary", we have observed a sharp crossover to the 
regime of predominantly good folding, in which the error probability is suppressed 
to a modest level ~ 10 — 20%. Interestingly, the behavior is nearly independent 
of the energy parameter a, which makes Rg/d the only major control parameter. 
Figure 5.3 shows the error probability as a function of this geometric parameter for 
tetrahedral clusters of different sizes. As this plot indicates, the misfolding behavior 
gets suppressed as Rg/d exceeds 1. 




Figure 5.3: Error probability P as a function of aspect ratio Rg/d for tetrahedral clusters with 
modified soft-potential and realistic DNA-mediated attraction. Each data point on the misfolding 
profile represents 100 trials. 
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As a further test of the robustness of the model, we consider a modified attractive 
potential which deviates from the linear Hooke's law. For larger forces we enter the 
Pincus regime [24], where the end-to-end extension of the polymer chain r is related 

3 

to the external tension / ~ ra. This tension law incorporates the excluded volume 
interaction of individual linker DNA with themselves, which was not previously 
considered. Remarkably, the major result for error suppression carries over, as 
illustrated in figure 5.4. 
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Figure 5.4: Error probability P as a function of aspect ratio Rg/d for tetrahedral clusters in the 
Pincus regime. Each data point on the misfolding profile represents 100 trials. 

Our description has several limitations which are to be addressed in future work. 
In particular, our discussion is only applicable to the limit of modest coverage of 
particles with buffer chains. This case of weakly overlapping adsorbed chains (known 
as the "mushroom regime") is drastically different from the high-coverage "polymer 
brush" behavior [21]. Nevertheless, our major conclusions appear to be rather 
robust. The condition for error suppression (i^^/c? > 1) is the same in both the 
harmonic and Pincus regimes. 

In the previous chapter we presented a detailed discussion of a method to 
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self- assemble DNA-coded clusters. Within the context of the beads and springs 
model presented here we can discuss the role that jamming plays in preventing the 
one anchor structures from assuming the minimal second moment configuration. 
We performed simulations of the assembly of optimal colloidal clusters up to n = 9 
particles by numerically integrating the particles' Langevin equations. As indicated 
in Figure 5.5, the hard sphere system gets trapped in a configuration with a larger 
M2 than the optimal cluster, whereas the soft-core system is able to fully relax. The 
jamming behavior is largely determined by the single control parameter with d 
the diameter of the hard sphere. Beyond the critical value ^ ^ .5 the jamming 
behavior is either completely eliminated, or greatly reduced in the case of larger 
clusters. 
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Figure 5.5: Plot of the dimensionless second moment ^ as function of time for n = 9 particles. 
Results are shown for the case of hard spheres and also for a system with a soft-core repulsion with 
geometric parameter ^ = .7. The dashed line is the theoretical moment for the triaugmented 
triangular prism, which is the minimal n — 9 structure. 
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5.4 Discussion 



Previous studies have demonstrated an experimental implementation of 
self-assembly in a DNA-coUoidal system [1], [26]. The approach in these studies 
differs from our vision of controllable self-assembly in two major ways: (1) because 
the "linker" DNA chains(~ 20nm) used are much shorter than the particle 
diameter(~ l/^m), the particles behave as sticky spheres, and (2) there is little 
diversity in particle type, where the structures result from interactions of a one 
or two component system. We would like to discuss issues related to a future 
experimental implementation similar to the modified Beads & Springs model. 

In our molecular dynamics simulations we assumed that (1) the desired group 
of particles has already been locahzed in a small region of space, and that (2) 
the interparticle connections have already been made. There are a number of 
experimental challenges associated with implementing the self-assembly proposal of 
our simulations. In another manuscript [61] we provide a detailed discussion of the 
localization problem, which is the first major experimental intermediary. 

After localization, the next step is to make the desired connections between 
particles within the cluster. To do so one can add short ssDNA with sequences saSb 
to link particles A and B. The DNA marker sequence sa for particle A is a sequence 
of nucleotides complementary to the sa portion of the linker sequence saSb- The 
hydrogen bonding of complementary nucleotides forms base pairs which join both 
marker strands to the linker, creating a DNA bridge between the two particles. 
After the interparticle links are formed, they should be made permanent by ligation. 
Since the spring constants of the above dsDNA chains are too small to drive the 
self-assembly of a desired cluster, we propose to melt them either by changing the 
temperature or pH. As a result, the dsDNA links will be turned into ssDNA with a 
much higher effective spring constant (due to the shorter persistence length) . This 
will trigger the self-assembly scenario similar to the one discussed within the Beads 
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& Springs model. Note that DNA entanglements may be effectively eliminated if 
the procedure is done in the presence of DNA Topoisomerases. 

5.5 Summary 

We presented a model of DNA-coUoidal self-assembly which exhibits a tunable 
jamming-unjamming transition. The combination of a soft-core repulsion with a 
type-dependent long range attraction provides a natural funncling of the energy 
landscape to the ground state configuration. This is to be contrasted with the case 
of protein folding, where under physiological conditions the interactions between 
amino acids are screened to several angstroms. Because this lengthscale is much 
shorter than the spatial extent of the native structure, large regions of the energy 
landscape are flat, which prohibits formation of the native state on the basis of 
tunneling alone. As a result the folding rate is necessarily limited by the diffusion of 
amino acid segments looking for their desired nearest neighbors [89], [90]. The fact 
that the potential in the DNA-coUoidal system is long ranged is essential, allowing 
us to avoid the pitfall of slow particle diffusion. 

Within the Beads and Springs Model, we obtained the jamming phase diagram 
for several modest sized tetrahedral clusters. We identifled a regime of parameter 
space with error-proof folding, and demonstrated the importance of introducing a 
soft-core repulsion. The original model was then adapted to include several features 
of reahstic DNA-mediated interactions. Although the jamming cannot be completely 
eliminated in the modifled system, we identifled a regime of predominantly good 
folding, and calculated the error probability for tetrahedral clusters. The jamming 
behavior is determined by a single geometric parameter Rg/d. We concluded 
by discussing prospects for an experimental implementation of our self-assembly 
scheme. 



95 



CHAPTER VI 



COOPERATIVITY BASED DRUG DELIVERY 

SYSTEMS 

6.1 Introduction 

Nanoparticle based drug delivery systems have attracted substantial attention 
for their potential apphcations in cancer treatment [7], [8], [9], [10], [91]. It is 
hoped that by selectively targeting cancer cells with chemotherapeutic agents one 
can reduce side effects and improve treatment outcomes relative to other drug 
delivery systems which do not discriminate between normal and cancerous cells. 
For example, many epithelial cancer cells are known to overexpress the folate 
receptor [92] , [6] , [93] , [94] , [95] . A nanoparticle with many fohc acid hgands will 
preferentially bind to cancerous cells. A recent study [4] of a potential drug delivery 
platform consisting of generation 5 PAMAM dendrimers with different numbers 
of folic acid found that multivalent interactions have a pronounced effect on the 
dissociation constant Kd- This enhancement is the signature for cooperativity of 
the binding, which should lead to a greater specificity to cancerous cells in vivo. 

In this chapter [96] we present a theoretical study of these key-locking 
nanodevices (see Fig. 6.1). We introduce the idea that there are kinetic limitations 
to cooperativity-based drug delivery systems. In vivo the finite timescale for 
endocytosis prevents arbitrarily high cooperativity in the drug delivery system. To 
begin we provide a detailed analysis of the in vitro experiments [4]. Although 
enhancement of the association is the signature of greater cooperativity, in this case 
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Figure 6.1: A picture of the dendrimer "key-lock" binding to the cell membrane surface. 

it is due mostly to non-specific binding of tfie dendrimers to tlie surface. Due to the 
finite time window of tlie experiments, only indirect support can be offered to the 
notion of enhanced cooperativity. In the latter half of the chapter we expand the 
notion of kinetically limited cooperativity to the system in vivo. The equilibrium 
coverage of nanodevices on the cells is related to the concentration of folate-binding 
proteins and the strength of the key-lock binding. We quantify the preferential 
adsorption of nanodevices to the cancerous cells, and discuss how kinetic effects 
prohibit arbitrarily high cooperativity in the drug delivery system. The implications 
of the work for designing new drug delivery vehicles with enhanced specificity to 
cancerous cells are discussed. 



We now consider a simple model of the nanodevice system. A dendrimer with 
a maximum of M keys (e.g. folic acids) interacts with locks (e.g. folate-binding 



6.2 Key-Lock Model 
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proteins) in the cell membrane surface. A simple order of magnitude estimate 
for M ~ 30 can be obtained from the ratio of the surface area of the dendrimer 
to the surface area of the folic acid. In this way we implicitly take into account 
the excluded volume effect between the keys. The free energy for the dendrimer 
connected to the surface by m key-lock bridges is [30] 



The dimensionless energy parameter A contains information about the binding 
energy of a single key-lock pair, and the entropy loss associated with localizing a 
dendrimer on the cell-membrane surface. An estimate of A ~ 17.5 can be obtained 
from the dissociation constant of free folic acid K^^^ using the equilibrium relation 
between the dissociation constant and the free energy change for the formation of a 
single key-lock bridge, K^^ — |^exp(— A). Here is the localization volume of an 
"unbound" key. Below we determine the value ^ ~ 0.2nm from analysis of the in 
vitro experiments, which was used to determine A. 

The measured association rate constant ka of the dendrimer with folic acid is a 
factor of 10'^ times greater than kji^ of free folic acid. Only a factor of m can be 
attributed to the dendrimer having many folic acids attached to it. Here m is the 
average number of keys attached to the dendrimer. This pronounced enhancement 
of ka is the primary evidence for non-specific attraction between the dendrimer and 
the surface. 



The non specific attraction sq accounts for the Van der Waals attraction to the 
surface and hydrophobic enhancement. The experimentally measured ka values are 
reproduced by a reasonable energy scale —eg ~ 7T (see Fig. 6.2). 

We provide a simple explanation for the experimentally observed dependence of 
the dissociation rate constant kd on m. The dissociation rate constant of free folic 



m 



—TmA. 



(6.1) 




(6.2) 
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acid k^j' ~ 10~^ [s~^] provides a characteristic departure time of l/k^^' ~ 30 hours 
for those dendrimers attached by a single key-lock bridge. Moreover, the departure 
time for multiple bridge states increases exponentially in A, for two bridges it is 
exp{A)/kl^^ ~ 10^ hours. Strictly speaking the relaxation is multiexponential, with 
time constants for each bridge number. However, the experimental kd values are 
well fit by a single exponential. On the timescale of the experiment, we will only see 
the departure of dendrimers attached by a single bridge. 

The experiment measures the departure rate of dendrimers which are connected 
to the surface by a single bridge, but are unable to form an additional connection. 
Consider a dendrimer attached to the surface by one key-lock bridge. If the 
dendrimer has a total of j keys, the probability that none of the remaining j — 1 keys 
can form bridges is (1 — ay~^. We now compute the probability a that a remaining 
key is available to form a bridge. In the vicinity of the surface the dendrimer is a 
disclike structure [97] with radius a ~ i.Snm. By rotation of the dendrimer about 
the first bridge, a key located at position p searches the annulus of area 27rp^ to find 
a lock. The probability of encountering a lock in this region is 2Tr p^ag, where the 
surface density of the locks cTo — lo^mr?' averaging over the key location we 
obtain the final result 

1 r 

a — — 2np^(Todp ~ ^aag- (6.3) 
a Jo 

Assuming that during dendrimer preparation the attachment of folic acid to the 
dendrimer is a Poisson process, the probability of a dendrimer having exactly j keys 
is Pjijn) = ex]){—m)m^/j\. The final result is obtained by averaging the probability 
that no additional bridges can form over this distribution. The factor of j counts 
the number of ways to make the first connection. 

oo 

kd = oo = 4°^ exp(-Q;m) (6.4) 

E J^j(^) 
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Figure 6.2: Top: Plot of the association rate constant (Eq. 6.2) ka[M ^] versus m. Bottom: 
Plot of the effective dissociation rate constant (Eq. 6.4) /cd[s^^] versus m. In the fit 

= 10~^[s~^] and a = 0.15. The experimental data points are taken from Figure 5 in [4]. 

The formula predicts an exponential decay of the effective dissociation rate 
constant with the average number of folic acids on the dendrimer, which allows for a 
quantitative comparison to the experiment (see Fig. 6.2). Using a ^ 0.15, we can 
determine the localization length ^ ~ 0.2nm for locks in the experiment from Eq. 
6.3. This estimate for ^ is physically reasonable, and comparable to the bond length 
of the terminal group on the dendrimer. 

Similar to the finite timescale of the experiments in vitro [4], in vivo the 
endocytosis time provides kinetic limitations to cooperative binding. In equilibrium 
the concentration of dendrimers on the cell surface n is related to the concentration 
of dendrimers in solution Cgoi through the association constant Ka = n/{aoCsoi)- 
Although it is tempting to use our in vitro results to define the association constant 
as Ka = ka/kd, this approach is only valid provided there is a single rate for both 
association and dissociation. Because the dendrimer can form multiple bridges. 



100 



there are many different rate constants. We present a partition function method 
which accounts for the multiple rate constants in the problem, and for the possibility 
that in vivo there is surface diffusion of locks. 

6.3 Calculation of the Association Constant Ka 



In this section we consider the reaction Uo + Csoi ^ n for which the association 
constant Ka is defined as 

Ka = (6.5) 

^o^^sol 

Here (To is the concentration of locks in the cell membrane, Csoi is the concentration 
of dendrimers in solution, and n is the concentration of dendrimers attached to the 
surface of the cell membrane. To proceed we construct a vector s of length M, 
which is a hst of the possible sites fohc acid can attach to the dendrimer. If a folic 
acid is present at site i we have Si = 1, and otherwise Sj = 0. The concentration 
of dendrimers on the cell surface n is proportional to the partition function of the 
system. 

m=l tjtjjt...jtp 



n — 



X(7(ri)---(j(r^)exp 



T 

Here a{r) is the surface density of locks on the cell membrane at position r, and 
A denotes the total area of the cell membrane. The energy £jj...p(ri, • • • ^Vm) that 
appears in the Boltzmann weight is the elastic energy penalty required to form 
multiple bridges. The point is that in solution the dendrimer is roughly spherical, 
but must flatten to a pancake like shape to form multiple connections with the cell 
surface [97]. 

The ensemble averaging is performed by assuming that during nanodevice 
preparation the attachment of folic acid to the dendrimer is a Poisson process. 
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In this case (sj) = ^ is given by the success probabihty that a fohc acid 
attaches to the dendrimer, and the m point correlator {siSj ■ ■ ■ Sp) = (;^) '"- In 
other words, the probability of attachment of a given folic acid to a terminal 
group on the dendrimer is unaffected by the presence of other folic acids 
up to an exclusion rule which has already been taken into account. If the 
interaction potential between locks in the cell membrane is V{ri^ ■ • ■ , r^) we have 
(cr(ri) • • ■a{rm)) = (cro)"^exp [— l/(ri^ • • • , rm)/T]. By performing the ensemble 
averaging we arrive at the result for the equilibrium coverage of dendrimers 
connected to the cell surface by m bridges. Since the lock which forms the first 
bridge can be anywhere on the cell membrane, without loss of generality we place 
this lock at ri = 0. The integrand is then independent of ri, and the first areal 
integration gives a factor of the cell surface area A. 

in) - "-^A Y: / '"^ HOMv,) . . .a{rj) x (6.7) 

m=l 



[n) = Csoi^ 2^ / ((7„) exp 



m=l 



m 

m 



-V{0,r2, ■■■ ,rm) 
T 



X (6.8) 



Y [^y^W |^^^_6o + s.,.,(0^r,,---,r^) 



l\ ^-eo^'^'' (ma,exp(A))- e^(-^) 1 
{n) = c,oiC exp ( — 1 ^^(^J^ X (6-9) 

^ ' m=l "5 



(V^(0, r2, • • • , Vm) + £:y...p(0, ra, • • • , r^)) 



T 

Here X^'' — ^I^d' — C''6xp(A) is the association constant of free folic acid which 
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has been measured experimentally. 
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(6.10) 



This gives the final result for the ensemble averaged concentration of dendrimers on 
the surface (n): 
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m=l 



n, 
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exp 



(6.11) 
(6.12) 



Here the elastic energy e^^^ is defined by: 



exp 



— e 



(m) 
eZ 



(i^r2 • • • (Pvj-^ 

t2(m-l) 



(6.13) 



X ^ exp 



{V{0 ,^2, - ■ ■ , I'm) + £^ij---p(0, 1'2, • ■ ■ ) I'm)) 

T 



Defined in this manner, exp(— eg^^/T) has a physical interpretation as the 
Boltzmann weight for the elastic energy of the optimal m bridge configuration. The 
membrane surface can only accommodate a finite number of locks in the vicinity 
where the dendrimer is attached [4] . As a result = for m > mmax since forming 
additional key-lock pairs would require deforming the dendrimer into configurations 
prohibited by elastic stress and steric hindrance. 

The calculation of the equilibrium coverage above is applicable with and without 
diffusion of locks in the cell membrane. In the regime of fast diffusion the locks are 
free to diffusively explore the surface. Their positions are ergodic variables, and the 
overall ensemble averaged equilibrium coverage counts the Boltzmann weights for 
different lock configurations. In the regime of slow diffusion, locks are immobilized 
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in the cell membrane. This is the relevant situation when the locks have phase 
separated into protein rich (lipid rafts) and protein poor phases. 



When kinetic effects are taken into account, these two regimes are drastically 
different. When locks are diffusing, the dendrimer is able to attain the maximum 
cooperativity mmax- After the dendrimer makes the first connection, it simply waits 
for locks to diffuse in the vicinity of available keys to make additional connections. 
In the absence of diffusion, the optimal configuration can only be obtained by 
multiple binding and unbinding events, the timescale for which is prohibitively 
long. This is the case for lipid rafts where the locks are immobilized similar to 
the experiments in vitro [4], and the dendrimer is unable to attain the maximum 
cooperativity. This is the kinetic origin of limited cooperativity in the drug delivery 
system. 

We now quantify the preferential attachment of nanodevices to the cancerous 
cells, taking into account kinetic effects. Let n„i denote the concentration of 
dendrimers attached to the cell by m bridges. We can construct a differential 
equation for rim by considering linear response to the deviation from thermal 
equilibrium n^. 



Here 7 is the rate for endocytosis [98]. The dissociation rate constant kj^ for 
breaking all m bridges is: 



The steady state concentration is the solution to = 0. As a result we obtain 



6.4 Kinetically Limited Cooperativity 




(6.14) 




(6.15) 
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Figure 6.3: The ratio of surface concentrations of dendrimers on cancerous to normal cells liili^ 



as a function of A with r = 10. The dotted line corresponds to an endocytosis time 1/7=1 [hr] 
and the sohd hue is I/7 = 10 [hr]. Here m = 15, m^ax = 4, ^ = 3[nTO], and CTq = 2 x 10^'^[nTO^^ 
^ir^ = ST for m > 3 bridges. 



the total coverage n of dendrimers on the cell surface in the following form: 



n: 



eq 



m=l 



m= 



1 1 + 



(6.16) 



We now have a means to discuss the preferential attachment of dendrimers to the 
cancerous cell. The folate binding proteins on the cancerous cell are overexpressed, 
i.e. if their concentration on the normal cell is cTq, their concentration on the 
cancerous cell is rao with r > 1. The value of r is determined by the biology, and 
cannot be changed by the experimenter. To quantify the preferential binding of 
the dendrimer to the cancerous cells we calculate the ratio of coverage on cancerous 
to normal cells "i'''^",-* . Values of this ratio greater than r indicate the nature of 

n((To) ^ 

cooperative dendrimer binding (see Fig. 6.3). 
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Figure 6.4: Single-stranded DNA (ssDNA) on the dcndrimer hybridize to the ssDNA attached to 
the fohc acid (FA) key. 

6.5 Design of an Improved Drug-Delivery Plat- 
form 

The current experimental scheme uses direct targeting with fohc acid 
{ApA — 17.5), which does not optimize the coverage on cancerous ceUs. By- 
decreasing A the drug dehvery can be tuned to the favorable regime. To do so, 
consider binding to the cell through an intermediary, perhaps single-stranded DNA 
(ssDNA). Instead of folic acid, attach many identical sequences of ssDNA to the 
dendrimer. Then, one also constructs a folic acid-ssDNA complex with the ssDNA 
sequence complementary to that of the ssDNA attached to the dendrimer. The folic 
acid will bind very strongly to the folic acid receptors on the cell membrane, leaving 
the unhybridized ssDNA as a receptor (see Fig. 6.4). Effectively one has replaced 
ApA with a new value Adna which can be tuned very precisely by controlling the 
length and sequence of the DNA. Due to the large degree of overexpression [92], 
this change substantially increases the ratio of dendrimers on cancerous to normal 
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cells. As indicated in Fig. 6.3, with r ~ 10 there is a 5 fold improvement over direct 
targeting with folic acid! 

6.6 Summary 

In this chapter we presented a theoretical study of a cell-specific, targeted drug 
dehvery system. A simple " key- lock" model was proposed to determine the effective 
dissociation rate and association rate constants of the dendrimers as a function of the 
average number of folic acids, which permits a direct comparison to the experimental 
results. The equilibrium coverage of dendrimers on the cell surface was calculated, 
and the differences between in vitro experiments and in vivo studies were discussed. 
The degree of cooperativity of the drug delivery system is kinetically limited. We 
quantified the notion of preferential selection of dendrimers to cancerous cells, 
and demonstrated that the selectivity can be enhanced by decreasing the strength 
of individual bonds. A particular implementation of this idea using ssDNA was 
discussed. 
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CHAPTER VII 



CONCLUSIONS AND PERSPECTIVES 

7.1 The Relationship Between CrystaUization and 
Jamming 

After the publication of [65] two recent studies reported the self-assembly of 
colloidal crystals using DNA mediated interactions [29], [28], [54]. Up until this 
point most of the studies reported aggregation of colloids into disordered, amorphous 
aggregates [1] or close packed crystals [26]. In these earlier studies polystyrene 
spheres of diameter d ~ were grafted with dsDNA of length L ~ 20nm. 
Because L < Ip c:i 50nm for dsDNA the DNA linkers behaved effectively as rigid 
rods. The essential change in the new studies was to work with flexible (L » Ip) 
DNA with contour length L > d. For this purpose Gang et al. [29] grafted g ^ 60 
ssDNA onto gold nanoparticles of diameter d ~ llnm. Interestingly, they found 
that as one varies the length L ~ Na of the DNA (a ~ A3nm for ssDNA [99]), 
the system goes from a disordered to crystalline configuration. Increasing the 
lengthscale for the repulsive interactions between particles provides a means to 
smooth the energy landscape to the equilibrium body-centered cubic configuration. 

It seems that an important control parameter in such systems is the ratio ^/d 
where ^ is the lengthscale for the repulsion and d is the hard sphere diameter. In 
fact it has been theoretically predicted that one can obtain different crystalline 
morphologies by varying this control parameter [55]. This control parameter is 
something we have already encountered in chapter 5. There ^ ~ i?^ was the radius 
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of gyration of the buffer DNA. Of course a direct analogy between the two systems 
is not possible, as the self-assembly of a modest number of particles into colloidal 
clusters is quite different than the crystallization problem. Motivated by their 
similarities, it is at least plausible that the explanation for which systems are able to 
attain their ground state configuration and which remain "jammed" in metastable 
configurations is qualitatively similar. 

To proceed with such a comparison we have to determine an appropriate 
definition of the control parameter. Here I propose that by defining ^ ~ one 
might be able to obtain a qualitative understanding of which systems crystallize and 
which remain disordered in terms of a critical value of the control parameter R/d. 
Here R is the so-called "coronal thickness" for a polymer brush attached to a small 
colloidal particle d <^ R. To determine it! we consider a simple scaling argument 
of Daoud and Cotton [100], [39]. By balancing the osmotic pressure against the 
polymer elasticity we obtain 

Ic'a'T (7.1) 

Assuming that the monomer concentration takes its average value c = {Ng)/V with 
V — (47ri?^)/3 yields the result for the coronal thickness 

R^g-^N^a. (7.2) 

We are now in a position to compare the value of the control parameter for the 
different systems studied experimentally [29]. 

Table 7.1: Control Parameters for Equilibrium Colloidal Crystallization 
~/V R/d Aggregate? 

18 0.48 amorphous 

30 0.66 amorphous 

40 0.78 amorphous 

50 0.89 crystalline 

65 0.95 crystalline 
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It is interest to note this approximation gives a critical value of the control 
parameter ^ > 0.9, which is very close to the critical value given for the self-assembly 
of colloidal clusters in chapter 5. The value = 40 is the mean for a system 
prepared with particles of type A for which Na — 50 and particles of type B for 
which Nb = 30. All other systems were prepared in a homogeneous fashion with 
Na — Nb — N. This analogy seems to indicate that the system with shorter ssDNA 
is trapped in a metastable configuration analogous to the jammed configuration 
for the colloidal clusters, and the introduction of a long-range repulsive potential 
provides a kinetically feasible pathway to the equilibrium configuration. One topic 
for future research would be a more detailed consideration of the jamming analogy 
in colloidal crystals beyond the crude scaling argument given here. Comparison of 
the results to the experiments of [29] is potentially problematic, since they increase 
the length of the DNA but keep the grafting density constant for all systems studied. 
As a result the properties of the polymer brush that forms on the particle surface 
likely vary for different linker lengths. 

7.2 Building Up Complexity 

Studies on DNA-mediated colloidal crystaUization have emphasized the 
importance of thermal annealing as a means to facilitate crystallization. Here 
we explore another idea which might assist in the equilibrium DNA-mediated 
self-assembly of colloidal crystals. The approach, which to date has received little 
attention in experimental studies, is to self-assemble crystalline structures in a 
hierarchical fashion (see Fig. 7.1). Consider a binary system of particles with two 
colors A and B. The DNA marker sequences are chosen so that AB bonds form 
as a result of DNA hybridization, but not AA or BB. To take advantage of this 
cohesive energy, the nearest neighbours in the resulting crystalline morphologies 
should be of different color [55]. These types of crystal structures are commonly 
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described as a lattice with a basis [101]. For example, one can consider placing 
A and B type particles at the points of a body-centered cubic lattice so that each 
particle has 8 nearest neighbours of the opposite type. This structure, sometimes 
called the cesium chloride structure, is obtained by taking a simple cubic lattice with 
a basis consisting of an A type particle at the origin and a B type particle at the 
center of the cube (a/2) (x + y + z). Here a is the lattice constant and the primitive 
lattice vectors are x=(i Q 0);y— (O 1 O), and z = ( Q 1 ) • Similarly 
one may obtain the sodium chloride structure where the particles are located at the 
vertices of a simple cubic lattice and each particle has 6 nearest neighbours of the 
opposite type. This crystal structure can be described as a face centered cubic 
lattice with a basis. 

In fact, the same is true of more exotic crystal structures, for example the 
diamond structure. This structure, the same as that adopted by gallium arsenide, 
can be described as a face-centered cubic lattice with a two particle basis. The 
diamond structure is of particular interest for its potential applications in the 
self-assembly of photonic band gap materials [74]. Since these crystal structures 
can be built up by as AB pairs, one could attempt to assemble these crystals where 
the fundamental units are not individual particles but dimers. The dimers could be 
assembled first and then used as the fundamental components. 

Alternatively one could consider grafting two different types of DNA onto the 
A and B type particles. DNA with sequence Ai is chosen complementary to DNA 
with sequence Bi and A2 is chosen complementary to B2. By varying the length 
of the recognition sequence, one can tune the melting temperatures Ti and T2 for 
the AiBi and ^42^2 complexes. In fact, the experimenter can control the relative 
grafting density of the strands on the particles [37]. By choosing a low grafting 
density for the A2 and B2 strands, the system will behave in an interesting fashion. 
For temperatures T > T2 the particles will be dispersed as monomers. As the 
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Figure 7.1: Building up a honeycomb structure using dimers. The honeycomb can be viewed as a 
two dimensional triangular lattice with a two particle basis. 

temperature is lowered Ti < T < T2 we expect to see the formation of predominantly 
dimers. Finally, for T < Ti we expect aggregation of large clusters, and under 
appropriate conditions crystallization. The hope is that this intermediate stage in 
the assembly might facilitate faster crystallization, i.e. provide a kinetically feasible 
pathway to the equilibrium ground state. 

This is just one example of how complexity can be incorporated into the 
DNA-colloidal system. Other examples have already been considered earlier in the 
thesis, for example the self-assembly of DNA coded clusters. There we saw the 
potential advantages of working with a multiple colored system, as opposed to most 
of the experimental studies to date which have worked with at most two colors. 
These ideas naturally lead us to the next topic, the inverse problem in self-assembly. 
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7.3 The Inverse Problem in Self- Assembly 

A natural extension of these ideas is to consider the so-called inverse problem 
in self-assembly. Presume that you have a particular target structure S in mind. 
For example, this structure might be a particular crystalline morphology, or clusters 
with prescribed geometrical features. Given S, the goal is to determine the set of 
components C who self-assemble into S. The DNA grafted colloids considered in this 
thesis, or combinations thereof, constitute a particularly promising component class 
due to the diversity and complexity of the potential interactions. This procedure 
requires knowledge of the interaction potential which governs the dynamics between 
components. Some progress has been made on the inverse problem for isotropic 
colloidal systems which interact through a pair potential [102], [103], [104], [105]. 
This work determined a class of designer potentials which give rise to a number 
of crystal structures, in both two and three dimensions. The drawback to this 
approach is that the designer potentials are generally quite complex. They contain 
a large number of numerical fitting constants, so it is not known how the designer 
potentials could be realized experimentally. 

A first step towards understanding the inverse problem in the DNA-coUoidal 
system is to determine the pair potential for two colloids grafted with DNA. Some 
steps in this direction include [26], [37], [106]. In general one would have to consider 
the effects of anisotropy. The major goal is to understand how the experimental 
variables (grafting density of DNA, DNA flexibility and length, hybridization free 
energy, etc.) affect the interaction potential. By tuning these variables one might 
be able to construct potentials which mimic the designer potentials determined 
through Monte Carlo simulations. Of course knowledge of the potential which gives 
rise to a particular equilibrium structure is only part of the problem. The free 
energy landscape is generally quite rugged, so one has to take care to provide a 
feasible pathway to reach the ground state configuration on experimental timescales. 
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Future studies which combine insights drawn from experiment, theory, and modehng 
will hopefully lead to a realization of the incredible technological potential of 
DNA-coUoidal systems. 
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APPENDIX A 



BRIDGING PROBABILITY FOR SCHEME B 



In this appendix we calculate Ce// for scenario B (see Fig. 2.1). Taking the z 
direction normal to the planar surfaces separated by a distance 2/i, an eigenfunction 
expansion [67] yields the following expression for the probability distribution 
function. Here /3 = -^^i ^-nd N is the normalization (unimportant for our 
purposes) . 

P(r,r') = NPx{x,:i^)PY{y,^)Pz{z,^) (A.l) 

P^ix^x') = exp(^-^(x-x')^) (A.2) 
Pyiy^y') = exp(^-^(t/-t/')') (A.3) 
Pz{z,z') = ^sm(^jsm( ^jexp[-/?n2] (A.4) 

n=l ^ ' 

To circumvent the Dirichlet boundary conditions, we start the chains at a small 
distance A from the planes. We then have r'^ = {x\, y\, A) and v'^ — {x'2, y'2,2h — X). 
We first expand the distribution functions to first order in A and then perform the 



116 



integration over all the possible hybridized conformations ri —r2. 



NttX 



Pi = P{ri, r[ = {x^, y^, A)) = -^Px{xi, x^)PY{yi, Vi) x 



(A.5) 



n=l 



P2 = P(r2,r'2= (x;,y;,2/i-A)) 



NttX 
2h 



Px{x2,x'2)PY{y2,y'2) X (A.6) 



^(-l)-msin exp [-/Jm^] 



m=l 



The integration over the xi and yi components are Gaussian, which leaves the 
integration in the direction normal to the plane. 



/ 



jy2^3^2^2 00 00 

PiP2S'{r, - r2)dh,dh2 = J^i-^rnmlnmih) exp [-P{n^ + m^)] 



12/i2 



n=l m=l 



(A.7) 



Inm{h) = J sin (^^^^) sin (^^^^^p") d^i = hSn,m (A.8) 
The second line follows since m and n are integers, using limj^o = 1- Defining 



Ar = — r2 the numerator in the expression for Ce// is: 



J P{ri,r[)P{r2,r'2)S\ri-r2)d'ridh2 

00 

^(-l)Vexp[-2/3n'] 



(A.9) 



12h 



3 / Ar 



4 R„ 



n=l 



We also make use of the following result. 

ANnXR'^. 



I 



(A.IO) 



k=l, k odd 
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This leads us to the overlap density for complementary, flexible linkers with x — 



Ceff 



c{x) = 



C[X) 



exp 



3 / Ar 



4 Ro 



-Stt 



£(-l)"n2exp[-2/3n2 



n=l 



64x 



(A.ll) 



(A.12) 



E exp 



When the chains are strongly compressed, ~ ^ » 1, the asymptotic behavior 
of c{x) is easily determined as the sums converge rapidly. The more interesting 
physical regime is as particles approach at separations large compared to the radius 
of gyration of the marker strands, which is /3 ^ 1. To extract the asymptotics in 
this regime we can massage the summation as follows. 

-1 d 



5^(-l)Vexp[-2/3n2] = ^^E(-l)"^^P[-2/5 

n=l 



(A.13) 



n=l 



2^exp[-8/5n2] - ^exp[-2/5 

. n=l n=l 

oo 

J2 exp [-(3k'] = J2exp[-(3k']-J2^xp[-4pk'] 



(A. 14) 



k=l, k odd k=l k=l 

Then note that for any even function of n, f{n) we have: 



oo ^ / oo \ 

n=l \n=— oo / 



(A.15) 



Finally we use an identity derived from properties of the theta function [107]. 



(A.16) 
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Expanding the sums using these identities, the result of some straightforward but 
tedious algebra gives the asymptotic behavior of the binding probability. 

2 



AGb ^ AGb + T 



AGb ^ AGb + T 



3 I A 
4 



3 / Ar 



4 V R. 



+ log (i?,^c.) + log (f) 
+31og(a;) + ^ 

2 

+ log (R^gCo) - log 

-2 log(a;) + 



9 / 3 
4 V T 



for X <^ 1 



for X » 1 
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APPENDIX B 



BRIDGING PROBABILITY FOR SCHEME C 



In this appendix we calculate Cg// for scenario C (see Fig. 2.1). If the method 
of images is used to construct the probability distribution function for the flexible 
linker DNA, as opposed to the cigcnfunction expansion used in Appendix A, we 
arrive at the following expression [67]. 



P(r, r') 



NPx{x,x')PY{y,y')Pz{z,z') 

-3 



(B.l) 



2Rl 



[z-z' - Anhf 





3 / /n2 




~ exp 




— exp 



— exp 



2RI 



[z + z' - Anhf 



2Rj 



{z + z') 



(B.2) 



Since = 3> 1 we only need the n = term in the expression for Pz{z, z'). 
Px{x,x') and PY{y,y') are the same as in Appendix A(equation A.l). Once again 
we start the chains at a small distance A from the planes. The majority of the 
hybridized conformations will have the planar surfaces separated by approximately 
the linker length L. To simplify the discussion we take Ar = from the 
beginning of the calculation, which corresponds to orientations of the rigid linker 
with a small component parallel to the surface. We then have r'^^ = (0, 0, A) and 
r2 = (0,0,L + A — A). We first expand the distribution functions to first order in 
A and then perform the integration over all the possible hybridized conformations 



120 



Pi = P(ri,r; = (0,0,A)) 



6XN 
E? 



Px{xi, 0)Py(|/i, 0)^1 exp 



2RI 



P(r2, 4 = (0, 0, L + A - A) = -—Pxix2, 0)Py(y2, 0) x 



P2 

■"-9 



(B.3) 
(B.4) 



(L + A - 2^2) exp 



-3(L + A - ;g2) 



To impose the delta function constraint we write: 



Z2 



y2 



X2 



5{\ri - ral - L)d?r2 



zi + L cos 9 
yi + L sin sin 
Xi + L sin 9 cos 
/ sin 9d9d(p 



(B.5) 
(B.6) 
(B.7) 
fB.8) 



The integrations over Xi and yi are Gaussian. There is no azimuthal dependence, 
so the 4> integration gives a factor of 27r. We also need the following integral: 

4A^7rAp2 



(B.9) 



We define e = ^ and z = ^. Since we consider DNA bridges with the rigid linker 
aligned with a small component parallel to the planar surfaces, the upper bound for 
the polar integration is given by 6* max — 7- 



Ce//(Ar' = 0) 



27C=^ f . , 



4 



Iz{9,eX) — J dz{l — cos9 + e — z)z 



h{9,e,C)d9 (B.IO) 
(B.ll) 





X exp 



~e{z' + {l-cos9 + e-zr} 



We first calculate the z integration 1^(9, e, (). Since we consider all of the possible 
conformations of the short linkers between planes, we can see that the upper bound 
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for the z integration is Zmax = 1 + e. However, since C ^ 1 the gaussian decay allows 
us to extend the integration to oo. We define f = z — ^ and the small parameter 
P = 1 — cos 9 + e. Completing the square gives: 



lz{0, e, C) ^ exp 



I - exp [-3CVT df (B.12) 



Since /3 is small we take the lower bound for the integration to 0. 



Iz{0,e,C) ^ exp 



1 f0- 



4C V 3 V 2 3C' 

We now expand the 9 integrand in terms of the small parameters 9 ^ e 

27 



:B.13) 



Ce//(Ar 



0) 



32V37ri?3 



exp 



CV_i 

2 3 



(B.14) 



/ 



X / ^exp 



d9 



9(^1-6-3) 
32V37ri?3 



exp 



4cv 

4 



3C^ 



(B.15) 



AGc( Ar' = 0) ~ AGc + T 



log (i?^c„)+ log ^^^J+ lev 
-log(.^-5f.) 

We are interested in the minimum free energy with respect to the separation between 
planar surfaces, determined by ^^^ci^r -o) _ q Performing the differentiation we 



(B.16) 



find e* = ^. 



AGclAr' = 0,e = e*) ~ AGc + Tlog 



( -2 

(62 — 

AGc + 4.24T + Tlog [L^RgCo] 



_ 9 
e 4 



(B.17) 
(B.18) 
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